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1.  Introduction . This  report  summarizes  briefly  the  research  carried 
out  under  the  Contract  N00G14-67-A-0094-003  1 on  the  calculation  of  super- 
cavitating  hydrofoils  near  a free  surface.  This  work  appears  in  six  pre- 
vious publications.  Two  of  these  were  supported  in  part  by  another 
Department  of  the  Navy  agency  (the  Office  of  Naval  Research*).  The 
present  report  concludes  this  series  of  publications  with  a users  list  of 
instructions  for  the  calculation  of  supercavitating  hydrofoils  of  finite  as- 
pect ratio  near  gravity-free  surface. 

2.  Publications. 

2.  1 "Nonlinear  calculation  of  arbitrarily  shaped  supercavitating 
hydrofoils  near  a free  surface"  by  O.  Furuya,  J.  FI.  Mech.,  1975,  Vol.  68, 
part  1,  pp.  21-40. 


2.2  "Exact  Supercavitating  Cascade  Theory"  by  O.  Furuya,  J.  FI. 
Engr.  (ASME)  1975,  pp.  419-429. 

2.  3 "A  Note  on  Three-Dimensional  Supercavitating  Hydrofoils", 
by  A.  J.  Acosta  and  O.  Furuya,  J.  Ship  Res.  , 1975,  Vol.  19,  No.  3, 
pp.  164-65. 

2.  4 "Three-dimensional  theory  on  supercavitating  hydrofoils  near  a 
surface",  J.  Fluid  Mech.  , 1975,  Vol.  71,  Part  2,  pp.  339-359. 


2.  5 "Numerical  Procedures  for  the  Solution  of  Two-Dimensional 
Supercavitating  Flows  near  a Free  Surface"  by  O.  Furuya,  1975, 
Conference  on  Numberical  Methods  in  Ship  Hydrodynamics,  Naval  Ship 
Research  and  Development  Center,  Bethesda,  Md.  , (with  partial  support 
by  Tetra-Tetra-Tech  Inc.  , Pasadena,  Ca.  ). 


2.  6 "An  Experimental  Study  of  a Superventilated  Finite  Aspect  Ratio 
Hydrofoil  Near  a Free  Surface",  by  O.  Furuya,  A.  J.  Acosta,  1976,  ONR 
Eleventh  Symposium  on  Naval  Hydrodynamics,  London. 

3.  Discussion.  The  basic  concept  of  the  present  work  as  outlined  in 
publication  2.  1 was  the  exploitation  of  the  non-linear  free  streamline 
theory  of  Wu  and  Wang  (see  e.g.  Wu,  1972  "Cavity  and  Wake  Flows"  Ann. 
Rev.  Fluid  Mech.  4,  243-284)  by  introducing  an  additional  method  of  de- 
termining the  constants  of  the  non-linear  equations  to  be  solved.  This 
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method  enabled  the  much  more  difficult  cascade  problem  to  be  tackled  with 
success  (publication  2.  2).  The  two-dimensional  non-linear  free  streamline 
calculations  are  embodied  in  a linear  finite  aspect- ratio  correction  in 
publication  2.4  for  flat  plate  fails.  Experiments  were  reported  in  publica- 
tion 2.  b which  show  that  the  present  theory  accounts  nicely  for  the  ob- 
served data. 

In  the  next  section,  the  computer  program  listing  together  with  a 
description  of  the  input  cards  for  the  finite  aspect  ratio  correction  effect 
(publication  2.4)  is  given.  The  important  geometrical  characteristics  of 
the  foil  are  given  in  data  cards  4 and  5 and  the  plan  form  shape  in  section 
4.  2.  In  the  present  program  listing  (section  4.  5)  the  chord  length  descrip- 
tion is  fixed.  The  plan  form  is  linearly  tapered  from  one  half  of  the 
semi-span  to  the  tip  where  the  tip  chord  length  is  one-half  that  of  the  base 
chord.  This  plan-form  distribution  may  be  changed  as  desired  following 
statement  number  1131  of  the  program  listing  section  4.  5.  An  example  of 
input  data  is  given  in  section  4.  3 and  the  output  from  this  data  is  given  in 
section  4.  4. 
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( ( $ 6 S ( S G G G 6 6 6 G fi  6 8 $ G E G £ 5 G 6 G C G C G 8 C 6 C G $ 6 G 6 6 5 G S G D $ 6 G u 6 o G G G G 5 G 6 G G 6 0 6 G « 5 G S G 6 G G S 6 b S 6 S ( G 
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| I i 8 I 4 4 8 ! I 3 8 9 3 5 ! ! f 3 8 8 i 8 8 3 5 8 3 3 £ 3 9 9 3 8 8 P.  9 a i)  it  b 4 R a 4 3 S 98  3 t i « £ M 3 S 3 8 H J 3 8 o 5 4 M S 8 J 8?  S 3 S J I 

9 9 9 9 5 9 9 9 3 3 


Card  No.  7 


NFSUR1  NFSUR2  MONK  NG01  NG02  NG03 

0 0 0 0 0 C C 3 0 0 0 0 G 0 3 C 0 0 0 0 0 C C Ij  0 0 0 0 0 0 o n 0 0 5 0 0 0 D 0 0 S U 0 0 0 0 0 0 0 0 0 0 P 0 0 0 G D 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 C 0 0 

1 I ) t 1 I ? i » it  ii  a u m ft  it  ii  it  »» n r.:« 11 » u u •>':  11:10  «•  <:  mi  <x; « *3 133551  ? ? r: ::  ;i  is  :>  :i  :s  ton  u ih<  cstct:  :o  n 77  n n » n 11  n »u 

11 1 it  1 1 n n 11 1 1 1 1 1 n 1 1 1 11 1 1 n 1 1 it  n n 11  M 1 11  n 1 11 1 1 1 11 1 1 nn  m n n 11  f 1 1 1 1 1 1 1 11  n 
22222222222222222222227222272222222222222222222272222722222222272222222222222222 

3 3 3 3 3 J 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 2 3 J 3 3 3 3 3 3.3  3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 J 3 3 3 3 3 3 3 3 

4 4 4 4 4 4 4 4 4 1 4 4 4 4 4 -14  4 4 4 4 4 4 4 4 4 1 4 4 4 4 4 4 4 4 14  4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 

555555555555555555555555:5555555555555555555555555555555555555555555555555555555 
6 6 6 6 S 8 6 6 5 6 5 6 6 6 6 6 6 6 S 6 6 G C S 6 5 5 C 6 6 6 6 E 6 t G 6 G 6 C C t 6 6 t*  £ 6 6 5 6 8 6 S E 6 6 G C S 6 6 E 6 6 C 6 S G 6 6 6 6 C $ B 6 8 6 6 S 
77777777777777777777777777777777777777777777777777777777777777777777777777777777 
II  8 8 8 8 6 8 J 8 3 5 8 8 « 8 f 8 8 8 8 8 8 8 8 8 3 3 8 8 3 8 8 S £ 8 3 SC  8 i 8 8 3 6 8 i t 8 8 8 8 8 3 S 8 3 8 8 C 8 8 8 I 8 3 8 8 8 8 I 8 8 J 8 8 8 8 S 8 


9 9 3 9 9 3 9 9 9 9 9 9 9 9 9 3 9 5 9 3 9 9 9 3 0 9 5 9 9 3 9 S 9 0 5 5 0 3 0 0 1 'J  9 0 3 9 3 9 9 9 S 3 9 5 D 5 9 9 3 9 3 3 5 9 9 9 9 9 9 0 9 9 3 9 9 9 9 9 9 9 

1 1 1 « j 1 i 1 lien  11  u 11  u i:  im>  is  Kin  tin  it  js  :i  .1 .1  ><  ; « n ii  >:  : " i.:i'o<i  'mu  « <5  «mi  < 1 53  st  si:  in  sis; vusti.cu.  ntu.  cii.'NU  10  n 11  a 11  nun 

0;.  I 


Card  No.  8 ~ No.  11  (1=1, -4) 

CIRCN(I) 

0 0 0 0 0 0 D 8 9 0 0 3 0 3 Q G C 0 0 C 0 0 0 0 9 0 C 0 0 1)  0 0 0 0 3 3 0 0 0 B 0 0 Q G 0 0 0 0 B 0 0 0 U 9 C 3 0 9 D 0 0 0 3 3 C G 3 G 0 0 0 D ? 0 0 CO  0 0 0 

1 J 1 i 5 1 ? 1 1 1:  M ’>  i>  m is  H w » 19  :ti  ii  :i  n j«  n :s  ;i  m m .1 ;/  n x ii  !t .i  n 10  « •; «.«  <>  hum  51  :j :<  ills  11  :i  suoti  o mu  <;  a »i  0 :o  n :i  :i  n it  u 11 mo 

u n n m 11  m u 1 1 \ 1 m m m 1 n 11 1 ! 11 1 11 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 m 1 1 1 1 1 1 1 1 1 11 1 1 1 11 1 

22222222222222272227222222222222222277222222222222222222222222222222222222222222 

3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 73  3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 

4 4 44  4 4 4 4 4 4 4 4 14  4 4 4 4 4 4 4 4 4 -14  4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 1 4 4 4 1 4 4 1 4 1 4 1 4 4 4 4 4 4 1 4 4 4 4 4 4 4 4 4 4 4 4 14  4 -I  4 
55555555555555555555505555555555555555555555555555555555595555555555555555555555 

6 6 S 6 6 8 G 6 8 8 6 6 G G S G G 6 8 6 G G 5 5 S 5 G 6 G G S S G G 0 G fi  G G 6 $ C G f,  G f.  5 6 6 C 6 G 3 5 G G b 6 5 0 6 G 5 6 6 S 6 b 6 3 6 G 6 6 G G 8 « S S 

7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 > 7 7 7 7 1 7 7 7 7 7 < 7 7 7 7 7 1 7 / 7 7 7 7 7 7 7 7 7 7 7 7 7 7 J 7 7 7 7 7 7 7 7 7 

8 8 8 8 8 8 8 8 3 8 8 3 8 8 3 8 £ 8 8 8 8 5 8 8 i 8 8 3 3 3 3 8 3 8 8 8 8 8 8 6 8 3 8 f.  8 8 G 8 8 8 8 8 8 8 8 fi  8 8 8 8 8 8 3 l 8 8 9 $ £ S 3 8 J S 3 3 8 8 I t 

9 9 9 9 9 9 9 9 9 3 9 9 C 9 9 9 9 9 9 9 9 5 3 5 9 9 9 5 9 9 9 9 9 9 9 G 9 9 9 9 1 9 9 3 9 G 9 5 ‘3  3 9 9 9 9 9 ? 9 9 9 9 9 5 9 ? 3 5 9 9 3 3 9 9 9 9 9 9 9 9 9 J 

1 1 1 t j 1 1 1 1 >»  >i  u n 1:  ;*  if  11  •’  i!  ?i  m?>  -i  7mj  ■.)  :*:<  -i ::  :• « • ;i  .imtiKiauuu  'i  si'.ii'i's'.'  'ifhv  ;•  * *1  ><  ■’  '•’•mu 


Card  No.  12~No.  15(1=1-4) 

PSSD(I) 

0 0 C 0 0 0 C 0 C 0 C 0 0 0 0 C U 0 C 0 9 0 G 0 0 0 0 n G 0 G 0 8 0 0 B 9 0 0 9 0 0 G 0 G C 9 0 9 0 0 0 0 0 C 0 0 0 0 0 C 0 8 9 0 0 0 0 0 0 0 C fl  0 0 0 0 0 0 0 

1 1 1 1 j < 1 1 1 it  n a 1)  m i'.  * n i!  11  »n  Hi) .’i » 11 1 1 >1  >0  ii  r 3:  1 j>n  17  31  i> o <1  <3  *1 11  n «. «t  4x1:11: 1: :i :■ 11  ;i  :m« i? 41  ti «s i{ r i«  n n 1 n .1  imimiii 

1 1 1 1 1 1 1 I I 1 11  I 1 1 1 1 1 11  1 1 I 111  I II  1 I 1 I I 1 I I 1 I I II  1 1 1 I I I 1 1 1 1 I 1 I I 1 1 1 I 1 1 I I 1 1 I 1 1 I 1 1 I 1 1 1 1 M 1 

2222722722272272 22222777272722272 77 ? 2 2 2 7 2 2 2 2222 2222222222222222 22222222222222222 

3 3 3 3 3 3 3 3 3 3 3 3 3 lb  3 3 3 3 3 3 J 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 

4 4 4 4 4 4 4 4 4 ' 4 1 4 4 4 4 4 4 4 1 4 <4  4 4 4 H 1 < < 4 4 4 4 4 4 4 4 4 M 4 4 1 4 4 8 4 4 4 4 4 4 4 4 4 4 4 4 4 4 M 4 4 4 4 4 4 4 8 4 4 4 H 4 4 4 

5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 S 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 

G i 6 5 6 6 6 G S 5 S 6 5 C G 5 S 6 5 f C 5 S S E 6 (■  t,  (■  G G 6 6 G 6 6 6 G 6 5 G 8 G G G 6 0 G 6 G G 6 S G 5 S G G S G G S G G 6 G 6 G 6 6 G S G 6 5 S S ( S 

7 7 7 7 7 7 7 7 1 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 / 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 1 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 

| 4 1 1 8 8 8 3 8 j 8 8 I 3 C 8 8 i 3 8 i I 1 1 S 8 4 3 ; i 2 J i a i j > 8 3 8 8 S 3 8 S 8 8 8 8 8 3 f 8 i 3 8 I 8 3 3 3 8 8 ? 8 3 8 8 88  8 1 8 I 8 3 3 G 8 I 

9 9 3 3 3 9 9 9 3 3 9 3 3 3 9 5 0 5 3 ? 7 5 9 *1  9 9 3 ? G 9 9 9 3 3 5 1 r 9 9 5 3 9 5 3 9 J 9 9 3 9 3 3 .>  3 9 5 3 3 3 9 ? 5 5 9 3 3 3 3 9 9 3 9 3 9 9 9 3 9 9 3 
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Card  No.  16 

SXSI(l)  SXSI(2)  SXSI(3)  SXSI(4) 

0 0 0 3 0 C C 0 0 0 0 0 0 0 3 0 0 0 0 fl  c C C 0 0 0 0 0 0 0 9 o 0 0 r.  3 0 0 0 0 0 0 C 0 0 5 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 3 0 0 0 0 0 0 0 0 0 0 0 U 0 0 
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44  4 41444144444444  4 444444  4 444 -I  4 '4  4 4444444444  1 444444  4 4 44444444  4 4444444444444444444 

9 5 5 5 S 5 5 5 5 5 5 5 5 9 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 C 5 5 5 S 5 3 5 i 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 S S 5 S 5 S 5 5 5 5 5 5 S 5 S 5 

$ 6 6 6 $ 6 S 6 6 S 6 6 $ f 6 6 6 6 6 C G b 5 G G G G > 8 G G 6 6 6 <•  2 8 6 G G 9 C G G £ G 6 0 i E G 9 2 b G S 6 9 6 6 6 C 9 C 9 6 6 6 R 0 6 C 6 6 6 6 6 G 6 C 

7 7 7 7 7 7 7 1 1 1 1 1 1 1 1 1 1 1 1 1 111  1 1 1 11  1 11  11 1 111  111  111  11  111  1 1 1 1 1 1 11  11 11 1 1111 1111 11111111  111 

8 8 8 8 8 8 8 8 8 1 8 8 S 8 9 8 ? 3 8 8 8 E S 8 S II  G 3 £ 0 3 8 3 t 8 3 2 5 3 5 8 8 5 8 * 8 8 f.  6 8 8 3 8 3 8 8 8 8 C 3 8 8 8 S 8 3 5 9 • 5 8 3 8 8 8 8 8 ( 8 8 


9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 5 5 2 9 9 3 9 9 9 9 9 3 9 9 9 V.  13  9 1 2 9 9 9 9 9 9 9 9 9 9 9 9 9 J 9 i 3 9 9 D 9 2 9 9 9 S 9 9 9 Vi  ? 9 9 9 9 9 9 9 3 9 9 3 

I I I « 5 6 7 • 0 If  II  II  O II  r.  16  'MJ  II  1)11  II I)  Ii  15  7.  II  717'  »•  7 la  55 K 1*  )-■  4141  II  '141 4'  .(  tl  17  4' 55  'I  51  s.  64  JO  56  S'  5.  51U  t>  W 51(1  IS  66  4 . 6)  1)  II  II  » 74  IS  M II  II  TO  u 


4.3.2  Example  of  Input  Data 


12960 


0.0360 

0(50010=0 


*10009? 


Ool 076667=  00 
0. 1037*5?=  00 
0.01  7n/.5  7= -01 

— ©itrMlttOC^-Ot 

0.25137^1)'  01 
!?.2'55e92=  oi 
0.29979*56  01 
0.2*279530  01 

0. 920 7 f'?ir-02~"'. 266037 3f  01  0.29690136  00  3.12777196-09 
-ri3?92^f,ft*^0?‘'0.2509/nT5-t>l— or246628  9F  -090i'23’8399?--<>4— 
0.B967‘ 31 =-02-0. 2631299=  01  U.2bb3139=  00  0 . 779 3092 =-09 
0.76196. ’9 =-02-0.255336 76  01  0.33522906  00  0.39200766-03 


1 These  are  the  data  obtained  from  the  case  of  a 
'Submergence  - 1 chord  and  cavitation  number- O' 


4-4  An  Output  Example  Corresponding  to  the  above  Input  Data  Set-ups 


AST  CONTROL  POINT  ON  SP4N  I S*****PA I /8»PA I * 3/64 

ALFA  ZFRC  = O.lJoOOOOE  02 HHH*  0.1000000E  01 

INITIAL  CIRCNI  11=  0.1076066E  00 
Initial  C I °CN I 21=  >).  10 j7 2536  00 

INITIAL  CIRCNI  31=  3.8170491E-0I 

INITIAL  CIRCNI  41  = <).62b8004E-01 
CAVIT.  NO  = 7.1289999E  00 
PSI/I  11=  0.2513700E  01 
P«IM  21=  0.2545892=  01 
P5IZI  3) = 0.2497925=  01 

P=tZt4)  * 0.2527938^  01 

m*N  wash  ANCLE  IN  OtCREfc=  0.2637746E  01 

POSITION  CF  SPAN 1 NO.  OF  TTER.«  l 

BIGS=  0.228571 4E  01  ISPAN=  l AR  = C.40009C0E  01 
PSIZ=  0.25 1 3700E  01 

ITERATION  NO. = 1 

— XIIT^  0.5297018E-02 

XI2) =-0.2660373=  01 
X I 31  = ".2464010=  00 

X I 4 ) * 0. 12  77719E-04 
X!l>=  C. 920  70l8  = -02 

XI  21=-0 .2660373=  01 

— X L?  1 ■*  0.2464010=  00 

X I 4 ) = 0. 1277719E-04 
Xi;>»  0. 932556 1F-02 
XI 21 =-0.2661556=  01 
X I 3 1 = 0.2470894=  00 
X 1 4 1 ■ 0. 394S363E-05 


ITERATION"  NOi*~t 


Xlll=  9.S325561F-02 
XI 21  = -0.26 61556  = 01 
X I 3 ) = 0.24  70894=  00 
X 1 4 1 * 0. 3949 363= -0  5 
X(ll=  0.9O39927E-02 
xt?i=i-'V29R4368c-or 
X I 31 = 0.2463559=  00 
X 1 4 I = 0. 72413500-05 


XIII-  0.5035927E-02 
XI21 =-0.29843680  01 
-Xf?;--  ''.2463599P  00 
X I 4 1 * 0. 7241350E-05 
XIII-  0.9193704=-02 
XI21  —0.2813456=  01 
X I 3 1 = 0.2467390=  00 
X 1 4 1 * 0.  74  4 10 1 7E— 0 5 


— ITERATION  NO.  « -3 


XUI*  0.  91  93704=-02 
XI ?> *-0. 281 3456=  01 
X I 3 1 » n. 2467390=  00 
X 1 4 ) * o. 744101 r = -05 
XII 1=  Co  9201959=-02 
X 1 2T*-*.  2809869*- 0T~ 
X I 3 I = 0.2467627F  00 
X 1 4 1 = 0. 7492128= -0  5 


XU  I » C»  92P  l959  = -02 
XI? 1 =-0.2809869=  01 
XIM-  rt.r4676?7F--00 
XI4)  = 0.  749"128=-05 
Xlll*  0. 920 3 1 73F-C2 
XI2)*-0.2RC83632  01 


xm«  0.24676616  00 
XI4I-  0.74869736-05 


ITERATION  NO.-  3 


0. 92031736-02 
■0.28083636  01  - 
0.24676616  00 
0. 74868736-05 
C. 9203069s -02 
■0.2808512'=  01 
o. 24676576  00 
O.  7480665*^05- 


ITERATION  NJi 


XIII*  C. 92030696-02 
X(2I — 0.28C8512P  01 

X I 7 1 - 0.2467657*=  00 
X (4 1 - 0.74886656-05 

XI  I J^—Cs 920 3 1105-0? 

XI  2) - -0.2800448'=  01 

xm-  0.24676586  00 
XI4I  = 0.  74864826-05 
Mil  »-9.  11920936-06 
*=121  * 0.4172325*= -06 

6-1  ? I • -<*.  5894  899c-94 

FI4I-  0.0 

P(  1.11-  0.  10326515  02 
P<  1,21  — <\  85543396-02 
PU, 31— 0.35482516-01 
PI1.4I-  0.0 

P I 2#  11-— >.  2371669^  02  

PI2, 2)  — 3.15237606-02 
PI  2, 3 60  385896-01 
PI2.4I-  *.  2 5786446  04 

PI  3. 11=  0.0  

PI3.2I-  **.91523306  30 

Pf3  r 37-— 0.  1 5969336-01 

PI  3, 4 1—0.  10000  30<=  01 
P 1 4»  1 1 — ?•  2 7275096  02 
PI  4,21—0.37352246-02 
PI  4, 31-  0.  12946  736  02 
PI4.4I-  0.0 

SX67M  1-— 0.9203 1306-03 

SXSII2I--C. 28084486  01 
SXSII3I-  0. 2467658E  00 
SXSTI4I-  0.74884636-05 

CIFCNI  II-  ). l>65345t  00 

CLIFT-  0.3281366L  00  CDF  AC*  0.87923656-01  CM0M6NT 

1 16^  0»  A6-C9F6PT5  -BAstO-ON-  Ol CI  O- -0.  37046615-00  ■< 

PPP-  0.24717616  00 
RCP*  0.24P41796  00 
Rro - i".  25050596  00 

pnp.  J. 2534  7076  00  

prp.  0.2*73553^  00 


DOWN  WASH  ANGIE  IN  DEGREE-  0.3323103E  01 

POSITION  C=  SPAN 

8IGS-  9.2285714E  01  ISPAN-  2 AR-  0.4003009E  01 
PSIZ-  9.25458925  31 

ITERATION  NO;*"t 

X  C 1 » = 9.82424885-02 

x 1 2» — 0.2593046=  oi 

X I 3 » » 3.2466283=  00 
X I 4 ) - 0.23203495-04 
XU)-  9.8242488=-02 

xt?)*-o.2303046rot 

X I 3 > =•  ?.?46*283E  00 
X t 4 ) - 0.  23263495  -04 
Xdl-  0.837 ?248=-02 
XI?) —9.  25039965  01 
X ( 3 1 - 0. 24 7 2 352  = 00 


0. 9 33  724  8E -02 
0.2508996=  01 
0.2472352=  00 
C. R3451 56? -C5 
O.8O77685=-02 
0.2752370=  01 
0. 2465568=  00 
0.  142  74  78F-04 


I TER  ATI  0*1  NG.«  2 


XI! 

XI2 

XI  3 
X I 4 

XII 

X 12 

X I 3 

-XI  4 Oi- 1409924*^-04 


0. 80776855-02 


XIII-  0. 9216053c-02 
XI  2)— 9.2640360=  01 

X I 3 ) = 0.2465304=  00 
X | 4) * 0. 1405524=-04 

~X  I t1r  0.  92 1 5 643® -*02 

XI  2) —0.264  7099'-  01 
XI 2) - 0. 2465245=  00 
X I 4 1 - 0.143  16523-04 


ITERATION  NO.-  4 


--  0.82156435-02 
0.^64  70999  Of 

- 0. 2465245E  00 
= 9. 143  1652  = — 04 
» 0. 82 1 8981=-02 
—0.2643114=  01 

- 0.2465335=  00 


XI4»'r-0.t424802E*0<r 


ITERATION  NO.-  5 


XID*  0.  32 189815-02 

X I 2)  —9.  26411 145  01 
XI ?) * 3. 2465335=  00 
XI4I-  >.14248025-04 
XU  I 0.  52  1 7994  = -02 

XI  21  — 0.  2644308=  01 
XI3J-  =.2465398=  00 
X I 4 ) - 0.  14265625-04 


ITERATION  NO.-  6 


0.821 79946-02 
■0.2644308*=  01 
0.2469308F  00 
0.  1426962F-04' 
0. 921 8 329*=-02 
0.2643873*=  01 
0. 2469317F  00 
0. 14262595-04 


ITERATION  NO 


XI  821  83 29<=-0 2 

XI21--0. 2643873=  01 
X I 3 1 - 0.2469317!:  00 
X ( 4) - O. 14262595-04 
XIII-  0.821 ei95c-02 
XI  2) —0.26440145  01 
X I 3I~= — 0.  24693145  OO- 
XI41-  0.  1426545*=-04 


ITERATION  NO.-  8 


XIII-  0.  82  1 8 195F-02 

X I 2) — 0. 2644014E  01 

XI  ■>)  - 0.  246S3146  00 
-X 141-  0* 14265456-04 — 

XU  1 = 0.821 82855-02 
XI  21  --0.  2643957*=  01 
X I 3 ) = 0. 24693165  00 
X<4)  = 0. 1426347E-04 
FID-  0.  47683726-06 
— f-<  2 1 — 0-.  47683->2*=-06 — 
FI31  * 0. 52809725-04 
<=141  = 0.9536743F -06 
PI!  , l)  = 0. 1099408F  02 
PI  1,2)  — \ 92477725-02 
PI  1*3)=" 7. 35394545-01 

-P(l»4t»-0.0 

PI2.1I-  U.2509952E  02 
P I 2. 2 I* -0.  1 74 5240*= -02 
PI?, 31— 0.48276306-01 
PI  ’.41=  ?.  1353662*=  04 
PI3.ll-  O.0 
-PF3v?> 

PI3.3) 

P I 3 , 4 ) 

PI4, 1) 

PI4.2) 

P 1 4 , 1 I 


*'.15456435  01 

■*'.  icnooooE  oi 
0. 2875328F  02 
'.41479545-02 
0.12938055  02 


SX5II11-  0.82132656-02 

SXSII2>»-0.26439j7E  01  

SX'1131-  0.24693166  00 

5X5II4I-  0.14263476-04  

riecni  21-  0.1325797fc  00 

Ft  I P T-~3T3i 475 72e-<W~«>**6--  -0*843  38445-0 FF*eW6NT-—flv54‘T74'>45—00-t-/f>--0¥  3732*648  01" 

Ltrr  ORAG  COcFFTS  3Ai60  ON  01 CLO-  C. 35536086  00  COO-  0.9521812F-01 

PFP*  J.  247  3 4226  00  

R*-P=  0.2485848F  00 

PrP-  C.  25*'4  743  6 00  

RTP-  '. 25364 l I F 00 

RFP--  0*  2475  2838-00 

PCP*  0.25239406  00 
AG1  II-  )• 1 I32367e  00 
AGI  21-  C*.  44189856-07 


AGI  31 « 0. 5153619E-02 
AGI  A )»  0. 1560A57=-07 
AGI  5).  0.ll82862?-02 
ACC  6)«  0» 1 695681E -07 

AG t ^T)«  0. 1659336?-02 

TGA.  0# 1 778718E  00  CL3D*  O. 3557A35E  00 

DOWN  WASH  ANGLE  IN  DEGREE  = O.310A72OE  01 

POSITION  OF  SPAN 

BIGS*  n. 181232&E  01  ISPAN*  3 AR*  C.SOOOOOOE  01 
PS  I 7*  0.2A97925S  01 

ITERATION  NO^*-t 

XI 1 I * 0. 8AA7928F-02 
X(2)*-0. 2631289s  01 
X ( 3 ) = 0.288313A?  00 

X I A ) * 0, 77A  3092s -0 A 
XU  | » 0.  8A6  7928F-02 

X I 21 • -'V  2671285s -01 

X|?)»  0. 28831 3A?  00 

X ( A I * 0. 77A3092F-OA 
XC1).  O. 85 8 1232E-02 
X 1 2) * -0. 26325SOF  01 

X I 3 1 * 0.2891010s  00 

XI  A)-*  0.  3176709'^OA  - — 


ITERATION  NO 


0.  85  8 1232F-02 
•3.  2632590F  01 
0. 2891010s  00 
C. 31 76709F-0A 
0. 832621 AF-02 
■0. 2923987s  01 
0. 2893216s  00 
0.50TAA32F-0A 


ITERATION  NO.*  2 


Oi  2922997  r-or 
'3.2883216?  00 
0. 507  A A32S-0A 
0.8A63100F-02 
% 2767331?  01 
■3.  288 738 IE  00 


0. RA63100? -02 

O.  2767331?  01 

P. 2CR7281F  00 
O.A87189l=-OA 


Xtn-*-r>.  8A9A9«9r 


0.2792653?  01 
0.28971A6?  00 
.l.*»983188?-0A 


XU)*  o.  8A5A569E-02 
XI ?T *•0.2783653?  Of 
X 1 3 I • 28  871 A6?  00 

XI*)-  O.A982188?-OA 
XIII*  rt.RA636A5?-02 
XI2I • -r. 2776263"  01 
X I 3 ) . ".288730A?  00 


-»t  A fa  -0.^8^8155' 


ITERATION  NC.«  5 


XU)  « 0.  8A606A5E-02 
XI2)*-". 2776263?  01 


0.2887304?  00 
0.4948155F-O4 
0. 8458689? -02 
•0.27786475  01 
O.288  7252«=-00- 
0.49611705-04 


iteration  no.*  6 


X C 1 1 * 0. P458689E-02 
. ' X I 2 I *-0.2778647?  01 
X(3I»  0.  2887252'=  00 
O — X 1 4h*-0.  496 1 1705-0*1 
xm»  0. 8459397=-02 
' X I 2 > 0.2777766?  01 
XI3|»  0.2887276=  00 
X (4 ) * 0. 495651 35-04 


ITERATION  N J.  * 7 


-*«-»  --Of  84593975-0-2 — 
XI 2) *-0.2777766=  01 
X t 3 » = 0.2887276?  00 
X 1 4 1 * O.4956513F-04 
XIII*  0.8459136=-02 
X(2l*-o. 27780875  01 
X 13  >*~o.  288  72675  00  — 
X ( 4 1 * 0. 49584635-04 


ITERATION  NO. * 8 


XIII  * 0. 84  59 1 365-02 
X(2»— 0.2778087=  01 
xm  * 0.2887267=  00 
— X !*->-*  O. 49584635— 04 — 
XIII*  o.pc,5<;248'=-02 
XI2l*-0. 2777954=  01 
X 1 3 ) * C. 2887270=  00 
X I 4 1 * 0.  49  577755-04 
FID*  0.59604645-07 

-5121— 0.-29802  32^-06 

F I 3 ) = 0. 120  5802=-03 
F|4)«  0.0 

PI  1,11=  0. 10803356  02 
P 1 1 * 2 1 *-0. 6675333E-02 
P|!t3)*-0. 3213135E-01 


0. 2318025E  02 
0. l70007l=-02 
o. 4 1 28795E-01 
0.38944705  03 
0.0 

-".00585455-00 

0.14570535  J1 
1000000=  01 
0.2202988E  02 
!?.  291  7921E-02 
O. 8918197E  01 


SX51I11*  0.8459248E-J2 

5X5II21*- o. 277 79 a4E  U1  

SXSII31*  0.2887270E  00 

SX51I41*  0. 49577755-04  

fI=CNI  31*  0.8086526c-01 

-CL  IF*-*  0.  3127329c-00-C08A6*— 0>-4374596F-Ol-C«4OH5N5— Or3759468e  50  tYi 

LIFT  0»AG  COEFFTS  PAiED  UN  Ul Clil*  0.35J3753E  00  Cl)0«  0.9460557E-01 

POP*  0. 2892C71 E 00  

PCP*  C.  290660  E 00 
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RrP*  0.29 31032E  00  

PCP*  0.2965722®  00 
ROP*  C.30H174E  00 
PTP-  0.3068066®  00 

AG<  II*  % H32367F  0© 

ACI  2)*  1.4418985E-07 

AG<  31*  0.5153619F-02  

AG<  4|«  0. 1560457F-07 
tG(  5 ) ■ 0.1182862®-02 
AGI  6).  ).  1655681F-07 

~AGI  J77»-  0.  155«5336e-02 — 

TGA=  0.1773718E  00  CL30*  U.3557435E  00 

OCNN  VA«H  ANGLE  IN  DEGREE*  0.3859832E  01 

POSITION  CF  SPAN 4 NO.  OF  ITER.*  I 

BIGS*  0.I468948F  01  ISPAN*  4 AR * C.4i)33900E  01 
P«I2*  0. 2527958®  01 

rr=RATioN  Ntnrr-i 

XIII*  0. 7614926E-02 
XI?) *-0. 25  53B66®  01 
X I 7 I * 0.335229C®  00 

X 14)  * •>.  7820076®-03 
X(1I*  0.7614S26E-02 

-XI  2) *-0.2553866®  01 

X 1 3 1 * 0.335229C®  00 
X ( 4 I * 0.3820076E-03 
X ( 1 I*  0. 771901 5E-02 

XI 2) *-3.2555338®  01  ' ‘ 

XI’I*  0.3360614®  00 

X*t  4‘t  *~  0.19312383-03 

ITERATION  NO.*  1 

XO  ) * 0. 7709C1 5E-02 

XI  2) *-0.2555338®  01 

XI 31*  0.3360614®  00  

X 1 4 1 * 0.  1931238E-03 

-X-tl>~  O.  74  515?7«—<V2 

XI?)*-''. 28 1 5305E  01 

XI3I * 0.3352327®  00  

X 1 4) * 9.2787148E-03 

ITERATION  NO.*  2 

XIII*  0. 74 8 1527E-02 

-X-m»“*.-2nt,7303»-<rt 

X(?I*  0.335’327®  00 

XI4) * 0. 278 7148F-03  

XIII*  0.7616018E-02 

X(2)*-0. 2674354®  01  " 

X I 3) * 0.  3357065®  00 

-rt4)-*-i>iP58*I1l3®-13 

ITERATION  NO.-  3 

XIII*  0.76"6018E-02  - 

X 1 2 I • -0. 26  74  354®  01 

X C 3)  * 0.3357065®  00  - 

X 1 4 ) • O. 2585 183® -03 

-Ttrr*-T»rr5»49it*=tj2 

XI2) *-0.2775381®  01 

X ( 3 1 * 0.3356342®  00  

X (4  I * 0. 26  8 1 287E-0 J 

ITERATION  NO.*  4 

XIII*  0.75845115-02 

-X 1 7 )'***.  270538 13“<7t 

X < 3 I * 0.3356342®  00 
X (41  * 0. 268I287F-03 
XUI*  0.  7596057F-02 


XU) 

X 1 2 ) 

X I 3) 
X (4) 

— xtrr 

XI2I 

xm 

X (4  ) 


1,2691542?  01 

O, 3356737?  00 

P. 2647419? -03 


ITERATION  NO.-  5 


*111-  1, 7596057C-O2 
X<  21  *-0,2691542'=  01 
X (3  )*  0, 3356737F  00 
X<4>  * 0. 264  74l9r-03 
XC11*  0. 7592C56P-02 
XI 2)  *-(>»  26966 11?  01 
X <31 3-0.-  3356  5 ?7*=  - PO 
X ( 4 ) * O. 266 1 108E-03 


X«  1»«  0.7597056? -02 
X<2)*-0.2696611F  01 
X < 3) * 0.3356577?  00 
•X  < 26<rl  lt)8«-0* 

X(lt*  0. 759371 7?-02 
X< 2) »-0. 2694553?  01 
X<?)«  0.3356643?  00 
X (4  I * 0. 2655464E-03 


I T5R4T 1 ON  NO. * 7 


X<  1 1-  0.  759-3717*-^2 — 
X<2»— 7.2694553?  01 
X«3»=  0.33566435  00 
X14>«  7. 26554645-03 
X<1 >»  0.75930285-02 
X ( 2 • *-9.2695400?  01 
X-<  31 *-  O.  3356613^00  - 
X ( 4 ) * 0.26577975-03 


ITERATION  NO.*  8 


0. 7593028E-02 
0. 269540Q*  01 
0. 3356613E  00 
-J.2657W5-03 
0. 7593293F-02 
3.26950475  01 
0.3356625?  00 
0.2656851F-03 


«!!*-?»  7593293*^-02 

X ( 2 ) *-0. 2695047?  01 
Xt’l*  0.3356625?  00 
X<4)»  0.26568515-03 
XU)*  75931855-02 
X<2)* -0.26951 90e  01 

-XT-*!---*.  3>5  6 679^-00 

XI4»«  0. 2657329F-03 
Ml)*  0. H93093e-06 
F < 2 ) * 0. 1370907F-05 
F<3)*-0.!289248e-03 
F<M*-0. 19073496-05 
*11  rH* — \ 1 1 420?5®— 0> 

P<1, 21  — 0.9C22161e-02 
PU  ,3)*-0.  29171905-J1 
P C 1 * 4 I » 0.0 
PI2.1I*  0. 2303421F  02 
P12.2>*-3. 1970594? -02 
-PI  **  SI*— V 3 5 464F-0 1- 

P < 2 . 4 1 » r. 7268225F  02 
P<x,il*  0.0 
P < 3 » 2 1 ■ 0. 5892462=  00 


■ 


p 1 3 * 3t > 3.1311190=  01  — - 

Pt3*4i*-i. 1000000=  oi 

Pt4,l»«-0.1835823=  02  

P « 4,  2 1 *-0.  2477034=  -02 

•P««.3t*  n«  63926346  Ot 

PI4,41»  0.0 

SX«ni)»  0.7593185E-02  

SXSII21*- 0.26951 90E  01 

SX5II31*  0.3356619E  00  

S X S 1 1 4 1 * 0.  26573<:9=-03 

1 *4 1 »~0.  621  74  16E-01 

CL I = T»  0.2955419c  00  CORAG*  9. 79 19984E-01  CM0M6NT*-0. 2297443E  00  L/D*  0.3732065E  01 

LIFT  DRAG  C05FFTS  BASED  ON  Ul --CLQ*  1.3336667E  00  COO*  0. 8940530E-01 

P6P=  0. 3362200E  00 

Rrp«  0. 3379090E  00  - - 

RrP*  % 3407494E  00 

— prp»-o.-344rt24-E-  00 

AGI  1»*  3. I 132367F  00 
AGI  2)*  0.44189S5F-07 
AG  I 31*  9.  5153619E-02 

AGI  41=  0.1  *69457= -07  

AGI  51=  O.t  182862F-02 

-AGI -61*  9i 1695681F-0T 

AGI  71*  9. 1659336F-02 
TGA=  9. 1 7787186  00  CL 30*  0.3557435E  00 
C1PC0U1-  0.13658816  00 

CIRC0I2I*  O.1026369F  00  - - 

CIRC0I31*  0.8999723F-91 

rtPCTH47-*-  0167209646-91 

PSSCI  11*  0.25109056  ul 

PS6ri  ?1»  0.25432605  01 

— PSSOI  31«  0.24955916 

P65PI  41=  0.25257986  ul  - 

02WN  WASH  ANGLE  IN  DEGREE*  0.2623338E  01 

— POMTICN  CP  -5-PAN-"”--  1— NO.— OF'  t T6RV=— t 

B IGS=  0.22357146  01  ISPAN*  l AR=  0.4000900e  01 

P6TZ=  0.25109056  01  

ITcbaTION  NO. = l 

X 1 1 1 = 0. 9203110e-02  -----  - - 

XI21«-0. 28984486  01 

— X I^r=-9. 2467658^—00 

X 1 4 1 * 9. 74884826-05 

XIII*  0.9203! 10F-02  - 

XI21—  0.2898448r  01 

X I 31  ■ 0.2467658=  00  

X 1 4 ) * 0.74884825-05 

— xti  r*-0i-9r?eo«7®-02 

XI  7}  *-''.  28953986  01 

XI31*  0.2466132=  00  - — - 

X I 4 1 » 9. 737I247F-05 


X 1 1 1 - 0.92260876-02 
-rt-?t*--9.-2p09398=- 
X I 3 1 * 0.24661326  00 

X 14 1 • 0.73712475-05 
XUl*  0. 922 3953=-92 
XI ?>*-9. 28119866  01 

XI  ’1  * ''.2466093=  00 
— X(41=-9. 73466T5P-05- 

I 

X I ! I ■ 0.92239536-02 
XI2l«-0. 2811986=  01 


ITERATION  NO.’ 


ITERATION  NJ.= 


X 1 3»  * 0.2466093E  00 
X * 0. 73466756-05 
XUI*  0.92232416-02 
X(2) *-0.28128365  01 
— X < 3) -0. 24660745  -«6 
X ( 4 1 * 0. 735 1 538F-05 


I TER  AT ION  NO. * 3 


X(ll*  0.0223241E-02 
X ( 2t  — 0. 28 1 2836*  01 
X(  31  - 0.2466074*  00 
-X-( 43*-Oy7 3515385—05 
X ( 1 1 * 0.022343 15-02 
X < 2 J — 0. 28 1 2638F  01 
X I 3 1 - O.2466C706  00 
X ( 4 ) » 0. 7349685E-05 


I T6RAT1 ON  NO. * 4 


XtM-*-Oo  9233431 F-02 

X I 2 J =-0.28126385  01 
> X ( 31 » 0. 24660786  00 
X ( 4 ) * 0.73496856-05 
X(ll*  0. 92 23383E-02 
X 1 2 1 * -0. 281 2709*  01 
- -Xf37*--O.2466O705-0O— 
X ( 4 1 * 3. 735000 2f -05 
Fill*  0.11 920935-06 
* (2) — 0. 1783139E-06 
F (3| >-0.65445905-04 
F (4) > 0. 1907349C-05 

P «-*  l 7—0^1^329496-0*-- 

°<l*2)>-0, 853831 7E -02 
P (1,31— 3.35498086-01 
P ( 1 # 4) * 0.0 
P<2,  1 >*  2371669E  02 

P ( 2 . 2 1—0.  1518257E-32 

P<  ?r-37~9.  483 39655-01- 

PI2.4I-  26274056  04 

P<3.1  )>  0.0 

P(3,2I*  9.91538946  00 
P ( 3 « 3 )=  0. 15987356  01 
P ( 3 * 4 1 — 3. t 0000006  01 


PI4, 21— 3.37297446-02 

PI4.31-  J.  12955036  02  * 

P<4,4)«  0.0 

SXSIllt-  0.92233636-02 

SXST (2)*-0.23127o9E  01 

SX  5 I ( 37— O«-2466lV78e-O0 

SX<!(4).  0.73500w2E-05 

C'FCNI  11=  0.  1065493d  00 - - 

Ct  IF Tr  0.32842201  00  CDRAG*  0. 8799994E-0 l CMOHENT— 0. 6237481 E 00  l/D«  0.3732070E  01 

tIFT  0PAG  COEFFTS  BASED  ON  U1 -CLO*  0.3707883*  00  CDO-  0.9935188E-0t—  

UPPEh.  CAVITV  SHAPE  LCb*R  CAVITY  SHAPE 

I«-J- 1-  XXU—OvO — — VYO-  0.-C 

ISJ*  3 XXU— 0.2522762E-03  YYU=  0. 7310498E-03 
!'J*  5 XXU— 6.  l36436«.E-0i  YVU*  rt.  181 1706E -02 

ISJ-  7 XXU>  0. 37257  7SE-0J  Y YU*  0.316  40616-02 
ISJ*  9 XXO»  0. 1279542E-U2  YYO«  0.4738891E-02 
ISJ*  11  XXU*  0.259501VE-02  YVU*  0.65100716-02 

353*-i  3 -XXIJ--0. 4*328096-02-  YY«—O.-04frl  Tl-5*-O-2 

ISJ*  15  XXU*  O.65O949UE-02  YVU*  0. 10583626-01 
ISJ*  17  XXU*  0.91 43 36UE-02  YVU*  0.12869176-01 
ISJ*  19  XXU*  0.  12255i)jE-01  YVU*  0.15314246-01 


XXL*  0.5929726E-01  YVL  —0.  1 1 192486-01 
XXL*  0.1472265E  00  YVL— 0. 28482726-01 
XXL*  0.23616756  00  YYL— (1.4341830E-01 
XXL*  0.3258732E  00  YVL— 0. 56982936-01 
XXL*  (1.4162390E  00  YYL— 0.6955123F-01 


XXL*  0.598 7 166E  00  YYL  — 9. 92412356-01 
XXL*  0.69074206  00  YYL — 0. 10291 2 16  00 
XXL*  0.  7837475*  00  YYL— 0. 1 1288 19E  00 
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ISJ*  21  XXII*  9.1586717E-01  YYU*  0. 17916526-01 
ISJ*  23  XXU*  0.2-37  ) V5oE-01  Y YU*  0.2767*176-01 
ISJ*  25  XXU*  0.26694626-01  YYU*  0.23590586-31 
ISJ*  27  XXU*  0.29966356-01  YYU*  0.26664185-01 

-I*J*  29  XXU*  0.35852845-01  YYU*  0.29838376-01 

ISJ*  31  XXU*  0.423893SE-U1  YYU*  0. 2 32 9 <4 eE -01 
ISJ*  33  XXU*  0.4961 494 E-0 1 YYU*  0. 368 6 2 72*-0 1 
ISJ*  35  XXU*  0. 57572256-01  YYU*  0.436022 36- H 
16 J * 37  XXU*  0.663383dE-01  YYU*  0.44520995-01 
ISJ*  39  XXU*  0.75375346-U1  YYU*  0.4,7621145-31 

7SJ*  41  XXU*_'0.36330596-Or— YYU*  "0.5292*^75-31 

ISJ*  43  XXU*  0. 9 77380j*-0i  YYU*  0.5  7429155-01 

TSJ*  45  XXU*  0.1131689E  00  YYU*  0.62146876-31 

ISJ*  47  XXU-  0.12370286  00  YYU*  0.67390935-71 

ISJ*  49  XXU*  0.13842916  00  YYU*  0.72275045-71 

ISJ*  51  XXU*  9.15444856  00  YYU*  0. 7 7 7 1 48 9E -0 l 

rrjs-^VKXU*  ‘0.171  37656“ 00  VW!”  0.8342 7976-01 

ISJ*  55  XXU*  0. 19784256  UU  YYU*  0.89434336-01 

rsj*  57  XXU*  0.2114965c  UU  YYU*  0.95^58685-31 

ISJ*  59  XXU*  0.2340  LI  46  00  YYU*  0.19242666  00 

ISJ*  61  XXU*  0.2585378E  UO  YYU*  0.19946976  33 

ISJ*  63  XXU-  0.28546146  00  YYU*  0.11692456  37 

ISJ*  -65"  XXU*  0.3149103c  00  Y YU* " O . 1 248 3396"01 

ISJ*  67  XXU*  0.34726686  00  YYU*  0.13324975  07 

I*J*  69  XXU*  0.3829321c  00  YYU*  0.14223175  OC 

ISJ*  71  XXU*  0.42239696  00  YYU*  0.15185476  07 

ISJ*  73  XXU*  0.4662696E  00  YYU*  0.16220726  70 

TSJ*  75  XXU*  O. 515318 UE  0 u YYU-  0.17341C56  00 

IS\7*  77  XXU*- P«5  7;)529sE"00  YYU*  0.  18567YCE"7<7-  — 

ISJ*  79  XXU*  0.6332-1616  UO  YYU*  0.19897486  70 

ISJ*  81  XXU*  0.7051080E  00  YYU*  0.21377735  J7 

ISJ*  83  XXU*  P0  788  705 je  Ou  YYU*  0.23735726  70 

ISJ*  85  XXU*  0.88761476  00  YYU*  0.24721216  07 

ISJ*  87  XXU*  0.10)74^2E  01  YYU*  0.2717966E  70 

trj-*“89  XXU* — hi  t OTirjE  U t- YYU*  0.  2772470*  "Ot 

ISJ*  91  XXU*  0.13555946  01  YYU*  0.32986716  00 

ISJ*  93  XXU*  9.1647946E  01  YYU*  0.3737414E  00 

Iej  = 95  XXU*  0.21 3961 jE  01  YYU-  P. 44297476  00 

ISJ*  97  XXU*  0.1553769E  02  YYU*  P. 1778749=  01 

ISJ-  99  XXU*  0.15893956  02  YYU-  0.11403855  01 

TSJ*Kn-XXO*-0il«597086"Ur  YYU*  0.  I753f736  "01 

UPPER  WAKE  SWAP6  IS 

XW(  II-  0.16597086  02  YK ( II*  0.17536736  01 

XW(  31*  P.15469995  02  Y»I  31*  0.90667976  00 

XW(  5) *-0.13643646-03  YW l 51*  0.18117066-02 

FP66  SUR FA  CP  SWAPF 

— RC^*-7i^ 4 tO\  T8t~00 

POP*  i.  24825886  00 
RCP*  0.2573456E  00 
PCP*  0.25330856  00 

-POP*  ".25715065  00  

POP-  7. 2627459E  00 


XXL*  0.87620556  00  YYL—0. 1223699E  00 
XXL*  7.  9695  51 4F  07  YYL— <7.  1 31 41 4 76  00 
XXL*  0.10633836  01  YYL--7. 14004726  00 
XXL*  0.11575606  01  YYL--0. 14829406  00 
"XXL*  7.12521056  OIYYL  — 7. 1561770F"00 
XXL*  0.13477706  01  VYL--0. 16371546  00 
XXL*  0.1442230*  01  YYL  — 0. 1 7092506  00 
XXL*  0. 153778  IE  01  YYL*-7.  1778203E  00 
XXL*  0.16336386  01  YYl—O. 18441 406  00 
XXL*  0.  17297906  01  YYL— 7. 1907167E  00 
" XXL-  "7. I82622'36—71  YYL-'-'li  19673806  00 — 
XXL*  0.19225296  01  YYL*-0.2024859E  00 
XXL*  0. 20198965  01  YYL=-0. 20756796  DO 
XXL*  0.21171136  01  YYL--0.21 3 1903E  00 
XXL*  7.2214573=  01  YYL  — 0.218 1583E  00 
XXL*  0.23122666  01  YYL--7. 2228768*  00 
"XXL  * "0.  241 0185'*- 01  YYL*— 0 .22734966  00 — 
XXL*  0.25083216  01  YYL=-0. 231 5795*  00 
XXL*  0.2676667*  01  YYL* -7. 23 5569 56  00 
XXL*  0.27052156  01  YYL--0. 23932066  00 
XXL*  0.28039606  01  YYL--0. 24283406  00 
XXL*  0.29728946  01  Y YL=-0 . 246  109 76  00 
-XXL  *"0.-3072ri2*-OT"  YYr^=0T249 T4716"  (TO 
XXL*  0.3171308c  01  YYL— 7. 25  19450E  00 
XXL*  0.3207776*  01  YYL--0. 25450076  00 
XXL*  0.3370411*  01  YYL* -0. 2568 l C 66  00 
XXL*  0.3477208=  01  YYL*-0. 25887056  00 
XXL*  0.35001615  01  YYL=-0. 2606741*  00 
" XXL*  "Oi ‘36702666-" OT  YYL—0.  26  22 17  86  0 O'  " 
XXL*  0.37705206  01  YYL* -0. 26347996  00 
XXL*  0.38005166  01  YYL  — 0.2644604E  00 
XXL*  0.3971450*  01  YYL  — 0.2651402*  00 
XXL*  0.40021186  01  YYL— 0. 2654997E  00 
XXL*  0.4102516*  71  YYL  — 0.  2655137*  07 
— XX  t*~Ov  4 2T38 366-^H"  YYL  * -7v?6-514^4f  "Yl  7 — 
XXL*  0.43048816  01  YYL  — 0. 264361 36  00 
XXL*  0.44760386  01  YYL*-0. 263085  IE  00 
XXL*  0.45.773736  Ol  YYL* -0. 26 1 22 1 56  00 
XXL*  0.46086656  01  YYL— 0.2586012E  00 
XXL*  0.47101016  01  VYL*-0. 25486516  00 
~ XXL"*"t7.  4 8tr298*-0 r~YYt*-Cri~24775746  00 — 


■ - • ‘-7  - ■ '! 


POSITION  OF  SPIN 2 NO.  OF  IT=R.«  t 

BIGS*  0.2285714E  01  ISPAN*  2 AR  = 0.4030000E  01 
PST 2*  3. 254326C5  01 

ITERATION  NO.*  1 

Xttl*  0.82182858-02 

X(2I*-0.2643957F  01 
XI3I*  0. 2469316=  00 
X<4)  = 0. 14263475-04 
X (t>  = 0.82182858-02 
X<2» *-0.2643957=  01 

24693165-00 

X<4»*  0.14263478-04 

XIII*  0. 92  89973F-02 - 

X( 2) =-0. 264=2518  01 
X(?)  = 0. 2468992=  00 
X(4|=  0. 1356554F-04 

1 T5RAT ION-NO."*  "1 

XU)*  0.82859735-02 
XI 2) *-0. 2645251F  01 
X 1 3)  * 9.  246 8992F  00 
X ( 4 ) = O. 13=65548-04 
X<1>*  0. 82 7 80216-02 

X <21’-0. 2658412*  01  -- 

X « 3 » = 0.  2469676F  00 

X ( 4)  = 0.  1374490=-04 

ITERATION  NO.*  2 

XIII*  0. 82  7 802 1F-02 

XI  21* -9. 26584125  01 

X 1-3 1-*  0.-24 6 86  765-00 

X( 4) * 0.  13  744908-04 
Xll)=  0.42783295-02 
XI2)*-0. 2654561=  01 

X ( 2 1 * 0.  2468774F  00 

X (4»  * 0.  1367738F-04 

1 T 5R  A 7 1 SN-«t/r--3 

X<1»=  0.52783295-02 
X<2)=-0. 2654561=  01 

X I 3) * 0. 2468774=  00 
X ( 4 | * 0.  136  77388-04 
XUI*  0.82771 145-02 

-Xt?»-— O*  7655984«-0-l 

X I 3 » - 0.2468743=  00 

X I 4 ) * 0.13701635-04  

ITERATION  N0.=  4 

X(l>*  0. 827 71 14F-02  

X(2)=-0, 2655984=  01 

-Xt-»»-*-0-i2468743c— 00 

X ( 4 ) * 0.  137C163F-04 

XIU*  827  7554=-02  - - - 

X ( 2)*  -0. 2655476F  01 

X 1 31  * 9.24687=4=  00  — 

X(4| * 0. 13693305-04 

1 f 58  A T MW  Ntr.  *— 5 

XIII*  3. 92775545-02 

XI  2 1 ■*-"*.  2655476=  01  

X I 3 » = 0.  2468754=  00 
X(4)=  9.  1 365  330=-04 

XI  1>*  •?.  82  7 7409t-02 

-X  I 2 T*-Oj  2655651*- 01 

XI 31  * 9. 2468749=  00 
X I 4 | * 0. I36560ft=-04 
Fill*  0. 596C464F-07 
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F(21*  0.53S4418F-06 
F(-»)*-0.l6l6478c-03 
F(4)  *-0.9536743c-06 
P<1»1)=  1094938E  02 

P(l»21=-r’.  9195752c-02"  

PI  l,31*-0. 35403905 -01 
P ( 1 f 4 I * 9.  0 
PI  2, 11*  3. 2500415F  02 
P < 2«  2 ) *~na  1728149E-92 
PI  2t  3 )*-'■>.  4828725F  -01 

PC?. 4)*-  0.  1410222F04 

PI3,1»=  9.0 
P ( 3 « 2 I = 0.91 49395E  00 
P ( 3 . 3 ) * 0. 1 549452c  01 
PI3,4  1*-r«. 1 1090905  01 
PI4.1  i = -9. 2 E6579 IE  02 

~P  I 4 , 2 ) *-9.  3 7709  l9F=tT2 

PI4,  3)*  •>.  1255468F  02 
P14.4I*  0.0 

SX SI ( 1 ) » 0.82774O96-U2 
«XSt(2)=-0.26556ilF  01 
SXS!  < 31  - 0.  ?4637*t9F  00 

5X91(41=  0.  1 J696U05-U4 

C|PCNi  2 ) = 0„  1027434c  00 

CL  1 FT  * 0.315591st  00  CDkAG*  0.8456212E-0 

LIFT  0«AG  COEFFTS  HAjED  ON  lit CLO*  .1.356 

UPPE8  CAVITY  SHAPF 


l C MCME  NT*-0« 5974817E  00  L/D*  0.3732065E  01 
302 7E  00  CDO*  0. 954706 1 F-01 

LOWER  CAV'TY  SHAPE 


t | f 

ISJ*- 

ISJ* 

ISJ* 

- 1 

*SJ  = 

1 

1 

[ ; 

! s J* 

1 * 

ISJ*- 

I [ 

ISJ* 

ISJ* 

1 I 

ISJ* 

1 XXU’  0.0  YYU*  0.0 

0-XXU=^0.209753Ut-OJ  YYU*  0. 6908395E -03 ~ 
5 XXU*-9.519782eE-04  YYU*  0 . 1 7 1 8 43  1 E-0 2 
7 XXU*  ".50  14182E-03  YYU*  0.  33 1 73 1 IE-02 
9 XXU*  0.145 3 409 E— 02  YYU*  0.4510004E-02 


XXL*  O.n  7YL*  0.0 

X X L“*~ 0.526 5 C 77E-TT1  Y Yt - -OV 1 6T943  8 E=OT 

XXL*  0.13391575  00  YYL*-o. 3227679E-01 
XXL*  0. 216955 lc  00  Y YL* -0 . 456 1 729E-0 1 
XXL*  0.2930516=  00  Y YL =-0. 5 7 7 27 34E-01 


13  XXU*  0.46)9 77 ^E -02  YYU*  0. 89 6 5 5 37E - 12 
15  XXU*- 0 i 63 34 56 56-02  YYU*  0.  I 1193 1 2P-=T«- 


XXL  * 0.3822165F  00  YYL*- 0. 689 504 7E-01 
XXL*  0.4660558E  00  Y YL*-0 . 7946545E-01 


"XXV*- 0^5504  S97F—O0  YYT- -9v8938235F"-ir0r" 
XXL*  0.6352499E  00  Y YL * - 1 .937 77 1 2e -01 


ISJ*  21  XXU*  0. 164995 3E-0 1 YYU’  9,17192146-01 
1'J*  23  XXU*  9®  70 6 2 36 ce -01  YYU*  0.  19 7 39?0C-0 1 
ISJ*  25  XXU’  O.2540  72  3E-O1  YYU*  0 . 2 3 52 6066-0 1 
ISJ*  27  XXU*  0. 3977557E-01  YYU*  0.25463935-91 
ISJ*  29  XXU*  0.  367666oE-0i  YYL.*  0.  295  5 5 1 7E -0  1 
1SJ*  31  XXU*  0.434l66uE-01  YYU*  0 . 3 1 8 1 2 58E-1 1 
ISJ*  31  XXU*  0.50764976-01  YYU*  0 . 752 104 5E -0 1 
I $ J=  35  XXU*  '■’«5fl3552«.E-0i  YYU*  0. 3973326E-9t 
T c J*  37  XXU*  0.67735UE-01  YYU*  0 . 4 2 5 2 1 75E-0 1 
1*J-  39  XXU*  o.  1745  755F-01YYU*  O. 4644 75 35 ’Ot- 
I c>  J * 41  XXU*  0.  8808  JOtE  -01  YYU*  0.  505032  l£-0  1 
!SJ=  43  XXU*  0. 99 6698 7 E-0 1 YYU*  0.548  52376-01 
! S J * 45  XXU*  0.112296VE  Ou  YYU*  0. 5935 5 26E -0 1 

ISJ*  47  XXU*  0. 126042  bE  00  YYU*  0.64.1  7285  = -0  l 

I'J*  49  XXU*  0.14 ''99846  00  YYU*  O. 69 1 1 8 1 7F -0 1 

T ' J*~  5 1 XXU* ~0. I S 72660E  OO'TYU*  0.  7420552*-(ir“ 
1SJ*  53  XXU*  O. 1 74961 VE  00  YYU*  0 . 796 5 1 8 3E-0 1 

1SJ*  55  XXU*  9. 1942100E  00  YYU*  0. 853 76 13E-0 l 

1SJ*  57  XXU*  0. 21 5 1 87bE  00  YYU*  0. 9 1 400 1 5E-0 1 

|sj*  59  XXU*  0.2390450k  00  YYU*  O. 9 7 7 49236-9 l 

! SJ*  61  XXU*  '■)•  2629935E  0 o YYU*  0.11445275  00 

153*-6  3 *XU*  0.7992 7366  Oil  YYU*  0.11154355  00 

t SJ*  65  XXU*  0. 32  )1 66iE  00  YYU*  0.1190650E  00 

ISJ*  67  XXU*  9.353018  »E  00  YYU*  0.1*706185  00 

TSJ*  69  XXU*  0. 3892089E  Go  YYU*  0.1155920E  00 


XXL*  0. 7204673F  00  YYL* -0 . 10 7 704 3F  00 
XXL*  0. 8060799E  00  YYL  *-•).  1 1620656  00 


XXL*  0.3920676E  00  YYL*-0. 124 31 706  00 
XXL*  0.97H4121F  00  YYL *-0.13206396  00 


XXL*  ff;  196^096?  Ot~ YYL*-CU1394T93F-00~ 
XXL*  0.  1 1 5? 105E  01  YYL *-0.14655246  00 
XXL*  1.  123P426F  01  YYL*-0. 1533298F  00 
XXL*  0.1327044F  )l  Y YL* - 1.  159 8 1 5 76  00 
XXL*  0.14 1 4947E  01  YYL*-0. 16602246  00 
XXL*  0.  1503  124F  91  YYL*-1. I 7 196  046  00 
- “X  X L-*-  Jv  159 1 5 6 5 F~  <71— Y Yt *-9yTT7 63  8 95  09 
XXL*  0.16802606  01  YYL *-U. 1 8306625  00 
XXL*  0.  1769198F  01  YYL*-0.  18824935  00 
XXL*  0.185837 3 E 01  YYL* -O.  19 3 19475  00 
XXL*  0.  19477745  01  YYL=-«1.  19790  74E  00 
XXL*  0.20373956  01  YYL*-0.20239?2E  00 


~XXT*-9V? T2 772 7 F 0 I~  YYL* -9. 29663? 6 F~00 

XXL*  9.2217263F  01  YYL=-U. 21069225  00 
XXL*  1.2307498E  01  YYL*-9. 21451346  00 
XXL*  9. 2 39  79?  3E  01  YYL *-0. 2 1 8 1 1 8 15  00 
XXL*  0.2438534F  01  YYL* -0. 22 1 50 74E  00 
XXL*  0. 25  79324E  01  Y YL* -0 . 22468 7 25  00 

X X L *'tT.?670  288F-91— YYL*-9v2276«rr9«^00 

XXL*  0.2761420E  01  Y YL* -0 .230 3862=  00 
XXL*  0.23527166  01  YYL*-0.2329131F  00 
XXL’  0.29441716  01  YYL*-C.2352205F  00 


E 


ISJ*  71  XXU*  0.529265u6  00  VVU«  0.1557250E  01 

ISJ*  73  XXU*  0.573791t,6  00  YYU*  0.1555556=  01 

ISJ*  75  XXU*  0. 5235707E  00  YYU*  0.16516016  00 

ISJ*  77  XXU*  9.5796035E  Oo  YYU*  0.17670356  00 

ISJ*  7«T  XXU*  0.653219*.E  00  YYU*  1.11935056  T)— 

ISJ*  81  XXU*  0.716176W6  UO  YYU*  0.2033356F  00 

ISJ*  83  XXU*  0.8010085E  00  YYU*  0.21897716  00 

ISJ*  85  XXU*  9,901 376 IE  00  YYU*  0.2367386E  00 

ISJ*  87  XXU*  0.10229656  01  YYU*  0.25731536  00 

ISJ*  89  XXU*  0.I1  75296E  01  YYU*  0.23153926  0(3 

ISj*~9l  XXU*-0.1»76209t  Crl— YYU*-  0.  31235156  09— 

ISJ*  93  XXU*  0. 166571 16  01  YYU*  0.3532099E  00 

ISJ*  95  XXU*  0.2171568E  Ol  YYU*  0.51715706  09 

ISJ*  97  XXU*  0.1578255E  02  YYU*  0.85369756  00 

ISJ*  99  XXU*  0.16155866  02  YYU*  0.85618236  90 

I S J = 101  XXU*  0. 1686392E  02  YYU*  0.861  )5376  19 

U88CR  KRKf~SHA 86  - 15 

XW I D*  0.16863925  02  YY>  I 1) 

XWt  31*  0.15720975  02  YWl  3) 

XWl  5 > *-0.51978285-05  YWt  51 

6866  SUB pAC 6 SHAPE 

BOP*  9.25728556  00 

25852786  00- 

RCP*  0.2506167E  00 
BTP*  1. 2535828E  00 
BCP*  (3. 2575692E  00 
BOp.  0.26233386  00 
All  II*  9.11207916  00 

*06  2 6* — (>«  536}  5996-07 

AGI  31*  9. 50592576-02 
AG(  51*  9.18615966-07 
AGI  51*  0.11332396-02 
AGI  61*  0. 19161966-07 
AGI  71*  9.15751966-92 

50»*—Vt-7605356r-3)0"Ct-31*-g.-35  29:0691 

OCWN  WASH  ANGLE  IN  OEGREt*  0.31665886  Gl 

POSITION  OF  SPAN- 

RIGS*  9.18123266  01  1SPAN*  3 AR*  0.50000006  01 
P S I 7=  0.25955916  01 

ITERATION  NO.*  1 

85552586-02 

XI2I=-9.2777955C  01 
X 13  ) * 0. 2887270=  00 
X I 5)  * 0.59  577756-05 
XIII*  9o8559258c-02 
XI2)*-0. 27779556  01 

XI  3)  * "0.-28  872-706-00- 

X I 5 1 * 0. 5957775c-05 

XI  1) * 9.1587396^-02 

X I 2 I * -'->.  2778960=  01 

X 1 3 I * 0.2885051=  00 

XI  5)*  0.58700556-05 


XXL*  0.3035779E  01  YYL=-0. 23730516  00 
XXL*  0.3127537=  01  YVL*-0. 23916226  00 
XXL*  0.32195506  91  YYL— 0.25078636  00 
XXL*  0.33115826  01  YVL*-0. 252 17006  00 
x x l*~  o.-3  a o 3 6 6 2 r - trr-  YYC* =0 .2  53  3tT<r6t-0Tr 
XXL*  (3.35959756  01  YYL*-0. 255 17886  00 
XXL*  0.35885156  01  YYL—0. 25577836  00 
XXL*  0. 3680S806  01  YYL*-0. 25508556  09 
XXL*  0.3773665=  01  YYL*-0. 25507676  00 
XXL*  0.  386656  7E  (31  Y VL*-<). 255 72 1 76  00 
X X L *-  Ov3959 38  26 -01-9fYt=“0 . 2*3978-78-00 
XXL-  0.50525036  01  YYL «-9 . 252 788 26  00 
XXL*  0.5155525=  01  YYL— 0.  25  105926  00 
XXL*  9.52387376  01  YVl*-0. 2386356=  00 
XXL*  0.5332021=  01  VYL*-0. 2351870=  00 
XXL*  0,55250876  01  YYL*-0. 228001 8E  00 


0.86105076  00 
0.72280326  00 
0.17185316-92 


I-T6*ATtGfrNtrv 


XI 1) = 0. 05873966-02 
Xl?I*-9. 27789606  01 
X 131  * 0.2886051=  90 
X 1 5 1 * 0.  58  700556-05 
xtl)*  0. 85  6 3719  - -02 
*121—1i-?7«»510*  01 
X I 3 ) * 0.28859376  00 
X I 5 1 * O. 58  73036c-O5 


* 9. 34817196-02 
=-0.2783510=  01 
« o. 2885037=  00 
= 0.4873076=-04 
= 0.,94339l3c-02 
—0.  2783218=  01 
= 9.28859398  00 
= 0.48  705  = 8 = -04 


Xlll=  0. 84839136-02 

Xf  ?) —0.2783218=  01 

X ( 3 » = 0. 2835939=  00 
X14)  = 0. 48  7D558  = -04 
Xll)*  O.  e483794  = -02 
X(2)— 9.27333828  01 
X I 3 > = 0.2=95934=  00 
Xt 4)*-  0.^871 104=-0V ~ 
f(ll=  0.17881 39c-06 
=12)= -0.238 4 186= -06 
Fill-'.  1496673=  -03 
F (4) =-0. ’9073496-05 
Pll  ,1)=  0.  1080 3 3 58  02 
P(l. 2t--P. 86548336-02 
P(l,3>— 0.32  142946-01 
PI  1,4)*  0.0 
P 1 2,  1 >*  ’a  2 3 1 83238  02 
P(2,2)=-n.l6927l3=-02 
P( 2, 3)=-0. 4130694= -01 
Pt2.4)*  3.3964761=  03' 

p 1 3 , l ) = n.o 

P (3,2)*  '’.9060509=  10 
P 1 3 , 3 ) * 0.  1458950F  01 
P 1 3 » 4 ) — P.  10030  00=  01 
P(4, l)*-0. 22029888  32 


31*  0.  8757072=  31 

4)»  0.0  " 

SX5J11)*  0. 84937*4 E-02 

SX«1(2)— 0.2783362E  01  - 

SXSI13)*  0.2985934E  00 

SXSM4  8=  0.49  711048=04 

-CIPfNI  3)*  0. 908913U-01 

CLIFT*  O.3130823E  00  CDRAG*  0.83889668-01  CM0M8 NT— 0. 372441 5E  00  L/0*  0.3732072=  01 

-LIFT  0 = 40  808  88 TF  BASED  UN  Ul CLO*  0.3534698=  00  COO*  0. 94 71 1 36F -01 

UPPER  CAVITY  SHAPE  LCHER  CAVITY  SHAPE  

ISJ*  1 XXII*  0.0  YYU*  0.0  XXL*  0.0  YYL*  0.0 

TSJ*  3 XXU'*-0.t493'7l'tE-O3~YYU*  P.'6766^rre-0T 
ISJ*  5 XXU*  9.  5*0 41  3 34E -04  Y YU*  0.  15594  346-02 
*«J*  7 XXU*  0 . 6564-296-03  Y YU*  0. 2 729 109=-02 

ISJ*  9 XXU*  0. 1645  16  7E -02  Y YU*  0.4  08  9 6668-32 
I 5 J*  11  XXU*  O.303777U8-02  Y YU*  o. 561  78056-02 
ISJ*  13  XXU*  0.4849 04 oE -02  YYU*  0.73912458-32 
I'M*-  15  XXU'*'“O.'909568 76-02~YYU*“  Oi9l 297078  *02’ 

18  J * 17  XXU*  O.979653uE-02  YVU*  0. 1 109  70  7E-01 
ISJ*  19  XXU*  0. 1 29724uE-01  YYU*  0.1319940E-31 
ISJ*  21  XXU*  9.t6S46WE-Ul  YYU*  0.15434336-01 
ISJ*  23  XXU*  0.?O94?7jE-Ul  YYU*  0. 1 T800 798-3 1 
ICJ*  25  XXU*  O. 255894  56 -0 1 YYU*  0. 2 0298 73E-0 1 
T'J*  2T  XXtr*-O.309160(jr-Pr  YYU’-n.22929096^OT 
ISJ*  29  XXU*  9. 369 551 4E-01  YYU*  0. 2569  329=- ) 1 
ISJ*  31  XXU’  9. 4344222c -01  YYU*  0. 2959403F-9 1 
l=J*  73  XXU*  0.507163OE-U1  YYU*  0. 3 16 3440E-01 


0.1175447=  00  YYL *-0.26803966-0 l 
0.1896594=  30  YYL— 0.  38405378-01 
0,2624310=  00  YYL* -0,489459  IE-01 
0. 33577888  00  YYL*-0.5871640=-01 
0.4096518=  00  YYL— 6786966E-01 


XXL’  0.5508284=  00  YYL— 0. 84671448-01 
XXL*  0.63407448  00  YYL  — 0. 92433516-01 
XXL*  0.7097272=  00  YYL— 0.9982181E-01 
XXL*  0.7857656=  00  YYL— 0. 1068664=  00 
XXL*  0.96217116  30  YYL  — 3. 11 35909E  00 
XXCT  0. 9399?69=~Ofr  YYL'r-0.  V20315 1 8 00" 
XXL*  0.10160158  01  YYL*-0. 126155  7E  00 
XXL*  0.1073423=  01  YYL»-0. 1320273=  00 
XXL*  0.1171137=  01  VVL*-3. 137642 JF  00 


ICJ*  35  XXU*  0.5872019E-01 
!SJ*  37  XXU*  0.675009SE-U1 
I6J*  ?<>  XXU*  9,771 1059E-01 
1SJ=  41  XXU*  0.87&0661E-01 

I5j*  43  xxu=  o.9905?8:>e-oi 

! S J*  45  XXU*  '•'.1115305E  Ou 
I=J*  47  XXU*  '»125189»;t  Ou 
ISJ*  49  XXU*  •1.1393479E  OU 
I S J = 51  XXU*  0.15S3980E  UU 
f S J . 53  XXU*  0.173353OE  00 

ISJ’  55-XXU=  <'.t9?345ie  00 

ISJ*  57  XXU*  )• 2133219E  OO 
1SJ*  59  XXU*  0.2355567E  00 
ISJ*  61  XXU*  9. 26.H499E  00 
1 5 J * 63  XXU*  0.2R70367E  00 
I6J*  65  XXU*  0. 3 164949E  00 

ISJ’  6 7 “XXU*  -9i  343856;»e  Co 

ISJ*  69  XXU*  0.3845214E  00 
!=J  = 71  XXU*  9.  V239301E  00 
ISJ*  73  XXU*  3.4678397E  00 
1 5 J * 75  XXU*  0.516066«E  00 
1 5 J*  77  XXU*  ''.5720472E  00 

1SJ--  79  XXU—9. 634679 jE  00 

1 5 J * 81  XXU*  0.7055209E  00 
1 5 J*  83  XXU*  79  l l’ObE  00 
TSj*  85  XXU*  0.8888355E  00 
!=J*  87  XXU*  0„10905lv£  ol 
JSJ*  89  XXU*  O.U53423E  01 

|S-J— 91-  XXU’  O.  135611  7c  01 

ISJ*  93  XX!J’  13.1640924=  01 
I ISJ.  05  XXU*  •"*.  2 1 3 34«;  jE  01 

I 5 J = 07  XXU*  0.  1 5 29345E  02 
ISJ’  99  XXU’  0, 1 56431 4fi  02 
ISJ*!01  XXU*  0. 1 631279E  02 

UPPC="  ViAkE"51tAP=— 1 S 

X'„(  11  * 3.  1631279E  02  YW  i 

Xw(  ?)*  0.1518906=  0*  YWI 
XWI  51*  0.  604 133  4E-04  Y„( 

F=F  = SUF  FACE  SHAPE 

RfP*  7.2890732F  00 

- ppp»  ?905  255E  00- 

PCP*  I. 2929674E  00 
POP*  1.  2=’64  348=  00 
onp,  0. 3009779E  00 
PTP*  ’066645E  00 
601  11*  >.  1120791E  00 

6*C"  2 1*  >.  4 3 fr 1 5995-97 

AC ( 31=  0. 5049247F-02 
aoi  4i*  >. ie6i596=-07 
AGI  61*  0.  1 133239E -02 
6CI  61’  >. 1916196c-07 
AGI  71*  0. 1 575106= -32 

*©#*  9.1  760535E-09  Ct30*-u.-3521u699-90 

D7WN  WA'H  ANGLE  IN  OEGKEc*  0.J822903E  01 

PUS l T 1 ON  C=  SPAN 

RIGS*  0.1468948E  01  ISPAfi*  4 AR=  0.4000300E  01 
PSI/  = 9,2525798=  01 

1 T=RAT10N  NO, * 1 

XI!  I*  9.  75931865-92 

XI ?l  *-0.2655190=  01 
X I 3 1 * 9. 3356619=  00 
X 1 4 1 * O. 76573295-03 


XXL*  0.1249145=  01  YYL*-9. 14301 12E  00 
XXL*  0.13’7434F  01  VYL«-0. 148 143 IE  00 
XXL*  0.14055935  01  YVL— 9.  1 530465E  00 
XXL*  0.  1434812E  01  YYL*-9. 1577284=  00 

XXL  * ' 04T 1 56 3 88  IT- 01  YTC*-v3.  1621950=-00 

XXL*  0.  1643191E  01  YYL=-9. 16645 16E  00 
XXL’  0.1722734=  01  YYL*-0. 1 705028=  00 
XXL’  0.  1802502E  01  YYL*-0. 1743528E  00 
XXL*  0.1882484=  01  YVL— 0.  1 780049E  00 
XXL*  0.  1962677E  01  YVL— O . 18 1 461  7E  00 

7(Xt-=-9i  294  3071E-01— rVt*-0.  l847252fl-00- 

XXL*  0.2123663=  01  YYL*-0.1877972E  00 
XXL*  0.2294438F  01  YVL— 0. 1906787E  00 
XXL*  J.2285399F  01  YYL»-0. 193 3700E  00 
XXL*  0.2366539=  01  YYL*-0. 19587 1 IE  00 
XXL*  0.2447848E  01  YVL  — O.  198  1 8 1 2=  00 

X X t * O i? 5 29  32  56  - 0 t-TYt'—  0^2002987^-00 

XXL*  0. 2619662E  01  > YL — 0 . 20222 l 3E  00 
XXL*  0.2692757E  01  VYL--O. 2039462=  00 
XXL*  1.2  7/4  703E  01  YYL—9.2054691E  00 
XXL*  0.2856796=  01  YVL— 0.206 7848E  00  ' 
XXL*  0.2939033F  01  V YL* -9. 20788 7 1=  00 

XXt*  0.302  34085-01— YYtr—t)v 208  76 769-00 

XXL*  3.  310351  RE  01  YVL— 0.  2094163=  00 
XXL*  0.3136559=  01  VYL— ".2098204E  00 
XXL*  0. 3269325E  01  YVL— 0.2099640E  00 
XXL*  0.3362215=  01  YVL — 0. 209826 IE  00 
XXL*  0. 34  35  22  3F  01  VYL— 0. 206  3791F  00 

XXt*'  Oi  1 8 344E  0 1- YYt—0v2065659e-99 

XXL*  0.3601574F  01  VYL  — 0. 2073929=  00 
XXL’  0.3684507=  01  V YL * -0. 205 7 1 86E  00 
XXL*  0.2766333F  01  VVl*-0. 20 342 33=  00 
XXL*  0.3851834=  01  VYL *-0.20021 10=  00 
XXL*  9.393514SE  01  VYL— 0. 1936526E  00 


0.4379 34 IE  00 
0.37196485  00 
0.1559404F.-12 


XIII*  0.75=3185=-02 
XI2)  — 0. 3695190?  01 
XI?)  * 0.3356619=  00 
XI4)*  O.265T329E-03 

XU)*  0. 7648766= -02  

XI  2)  —0.  269  6233=  01 
X I 3 ) * 0,3356259=  00 
X|4)  * 0#  2569319r-03 

ITERATION  NO.*  1 

xm*  0.76487669-02 

~ XI7)  — Oi  269  6 233*"  Ot 

X I 3)  * 3356259=  00 

X I A) * 0. 256S319P-03 
XU)*  0. 7638961F-02 
XI2)*-0. 2708599=  01 
X | 3 ) * 0.1355879=  00 

— X(A)*  <?~  75956019-03 

ITERATION  NO.*  2 

XII) * 0. 7638961E-02 
X I 2 I *-o. 27085999  01 
X 1 3 ) * 0.335*879=  00 
XI4).  o. 25996015-03 

TCfl)*  0.  7642299=^02  ~ 

XC2) *-0.2704488=  01 
X I 3) * •>„  3356017=  00 
X 1 4 ) * 0. 2584558  = -03 

" ITERATION  NO.*  3 

XID*  0.76422999-02 

— XI?)  *-0.2704488'=-01 

X ( 3 1 * 0.  3356017=  00 

X I 4 I * 0. 25345589-03 
XID*  o„  7640913E-02 
O XI2)*-0. 77061929  01 
X < 3 ) * 0.3355960=  00 

- -Xt4t*  0r25  99O53r-O3 

ITERATION  NO.*  4 

XU)=  9. 764C913E-02 
XI2) —0.2706192=  01 
X 1 3 ) * 0.335*9605  00 
X 1 4 ) * 0. 2589053E-03 

— Till T*-0i 764 1423=-02 

X I 7 ) *-'>.2705498=  01 

X 1 3) * 0.3355978=  00 
X 1 4 1 = 0. 2587487E-03 

ITERATION  NO.*  5 

xm*  0.76414239-02 

xm-*--- 0.2705498=  Ot 

XI?)*  0.335=978=  00 
X 1 4 ) * 0. 25  8 7487=-03 
Xll)=  '>.7641241=-02 
■*  Xl?)=-9.  3705738=  01 

X I 3 1 * 0.3355973=  00 

--  Xt4)*-0.?5Rer>22«=T» 

|T  = RATIC1N  NO.*  6 

xm*  0.76412419-02 
XI2)»-0. 270=738=  01 
XI?) * 0.3355973=  00 
X 1 4) * 0.2588022=-03 

— XI  It  76413799-02 

XI2)*-0. 7705656=  01 

Xl>)  * 0. 3355976E  00 
X 1 4 ) * 0.25875779-03 


Fill*  O.3344650r-06 
FC2l*-0. 2980232=-<36 
F(3»*  0, 733 13  71C-04 
FI4)  *-0.953674 3= -06 

PIl.lI*  -3.  11372576  02 

P11.2 J*-0.8976031F-02 
PC1.3»*-9.29l7764E-91 
PU.41*  0.0 
P t ?t 1 I*  ’'.2293587E  02 
P ( 2 • 2 1 *-'>.  1953279=  -02 

P12,-3r*-«.355?1548-Ot 

P * ? . A » = 1.  7461526F  02 
PI3,1)>  0.0 
P C 3 1 2 I « 9.88965465  30 
P ( 3«  3 1 * 0.  1314485F  01 
PI  3.4 )*-0. 10000006  01 

-PI4,t-»*-9i  18310556-52^ 

PI4»?)«-o,.  229101 5E -02 
P<4,3>  = 0.  6251791F  01 
P ( 4 , 4 ) * 0.0 

SX?[|1I»  0. 7641379E-02 
$X$I(2l»-0. 27056566  01 

5X51(3) * -0. 3 355976  E 00 

SX5K4)*  0.2587577E-03 

C I F CN ( 41*  0.622S5nit-01 

CL  [FT*  0.2962440fc  L’O  CORAG-  9.79378015-91 

LIFT  09 AG  C06FFTS  BASED  ON  U1 CLD*  0.3345 

UPPEh  CAVITY  SHAPE 

IS** fr  XXUx-Oi'V--  VYU*  0.0 

I9J=  3 XXU*-0.680935dE-04  YYU’  0.54370226-93 
ISJ*  5 XXU*  0.22017666-03  YYU=  0.1 3579656-92 
ISJ*  7 XXU*  0.3945127L -03  YYU*  0. 23739515-02 
I=J*  9 XXU*  0. 19629 1 IE -02  YYU1  0. 39671  755-92 
ISJ*  11  XXU*  0.  343794(jE  -02  YYU*  0. 490 1 5 4 3 = -'12 

> +S**~1"3  XXU»-Oi933491  /E-02 - YYU*  0.  6 3690695-9-2- 

I = J*  15  XXU*  0.767 1343E-u2  YYU*  O. 706 1 2 55=-02 
ISJ*  17  XXU’  0. 10466696 -01  YYU*  0. 94 72 3 7 7E -02 
ISJ*  19  XXU*  0. 1 3742356-0 i YYU*  0.11498656-01 
ISJ*  21  XXU*  0.17521756-01  YYU*  0.13437695-01 
ISJ*  23  XXU*  0.21830456-01  YYu*  0. 1 548e2lE-0l 

15-J-*  25-  XXU*-i>s2669631E-Ol  YYU*  0.  17649836-91 

ISJ*  27  XXU*  0.3214969 1-0.  YYlJ*  0. 1992295E-0 1 
!SJ»  29  XXU'  0. 38223766-0 1 YYU*  0.22308646-01 
ISJ*  31  XXU*  0.4495472E-C1  YYlj*  0.24008626-01 
ISJ*  33  XXU*  0.52332296-01  YYU*  0. 27425236-01 
TSJ*  75  XXU*  0.605500ufc-Ui  YYU*  0.30161436-01 

I-5-J*-  37-XXU’-0.-695O5795'Ut  -YYU*-  <V3  3i>?07  ?e— 

ISJ*  39  XXU*  0. 797073  76 -01  YYU*  0.36007475-01 
ISJ*  41  XXU*  0. 89993306 -01  YYU*  0. 39 1 2 ( 44E-0 1 
ISJ*  43  XXU*  0. 1 01 658 oE  0 u YYU*  0. 4 2 3 8 3 3 1 E-0 1 

ISJ*  45  XXU*  0.11435546  Ou  YYU*  0. 4578454E-01 

ISJ*  47  XXU*  0.12U1695E  OU  YYU*  0.40337516-01 

I-5-J»-49  XXU-*-  0.14319196  Ou  YYU*-0. 530505 76-93- 

ISJ*  51  XXU*  0. 1 59525oF  Gu  YYU*  0.56933296-01 

ISJ*  53  XXU*  9.17728465  Oil  YYU*  0. 6099652= -J1 

ISJ*  55  XXU*  0.196603VL  UU  YYU*  0. 652 5266= -0 1 

ISJ*  57  XXU*  0.21763406  00  YYU’  0.6071*986-01 

ISJ*  59  XXU*  <)•  240550 aL  tu  YYU*  0.74402816-01 

|0-J-*-6 1"  XXU*  0.2655 56 56  UV  YYU*  0.  197321  l«— 9t~ 

ISJ’  63  XXU*  0. 29289 1 36  uO  YYU*  0.84525946-01 

!SJ*  65  XXU*  O. 322837U6  UU  YYU*  0.50010056-01 

ISJ*  67  XXU*  0. 355730,»E  GU  YYU*  0.9  5814946-01 


XXL*  0. 3796023F-01  YYL*-0. 1527474F-01 
XXL*  0.95336686-01  YYL=-0. 25 589 196-01 
XXL*  0.1534C446  00  YYL*-G.3448987F-01 
XXL*  0.2123214E  00  YYL*-’>. 42577036-01 
XXL*  J. 27112646  00  YYL=-0. 50076756-01 
— XXL-*-Yhr3306  8a2F-00-YTt=— >ii57106-t-7f-0t- 
XX L - 0.39065325  00  YYL=-0. 63737696-01 
XXL*  0.45102105  00  YYL*-0. 70021 155-01 
XXL*  0.511 7624E  00  VYL*-0. 75992296-01 
XXL*  0.57285905  00  YYL*-0. 81678876-01 
XXL*  0. 6342940=  00  YVL*-0.8710295<=-01 
— XX3r*— 3 s-6960521-^-00" -VYt*~0  r»22#2  2-46-  Ot- 
XXL*  0.7581192=  00  YVL*-0.9723157E-01 
XXL*  0.8204823F  00  YVL*-0. 10 196365  00 
XXL*  0.8831293=  OJ  YYL*- 3. 10648865  00 
XXL*  0.94604905  00  Y YL= -0. 1 1081 59=  00 
XXL*  0,  1009 23 IF  01  YYL *-0.11496355  00 
— XYt~*  9yt'>  72t<T5P  *Yt  TYt*--hi  l 18508 18--0» 
XXL*  0.1136341=  01  YYL*-0. 12268546  00 
XXL*  0.12002526  01  YYL*-0. 12629086  00 
XXL*  0.12643885  01  YYL*-0. 12972886  00 
XXL*  0.  13287426  01  YYL  — 0. 1330030=  00 
XXL*  0.  1393307=  01  YYL«-().  136 1 1706  00 
— XXlr*-«)T 1-4580  Tfr*-0t-YYl*-0.t  350736*-  Of7 
XXL*  0.15230416  01  YVL*-0. 141 8747F  00 
XXL*  0.15831975  01  YYl*-0. 14452265  00 
XXL*  0.1653537=  01  YYl—0. 14701875  00 
XXL*  0.1719055=  01  YYL«-0. 14936355  00 
XXL*  3.17847476  01  Y YL» -0. 15 l 5580=  00 
— XX  t*-5  .1 8503OT6  01  YYt-*~0yt336O2i-6-(W 
XXL*  0.1916630=  01  YVL*-0. 1 5549525  00 
XXL*  0.19828106  01  VYL--0. 15723676  00 
XXL*  0.20491446  01  VYL— 0.1588251E  00 


IGA*  9.1778/186  00  CL30*  G.3537435E  00 


29 


ISJ*  69  XXU*  0.39197816 


ISJ. 

ISJ* 


71  XXU*  0.632977.sE 
73  XXU*  0.676665iE 


ISJ.  75  XXU*  0.526659*.E 
T'J*  77  XXU*  0.58252C5F 
ISJ*  79  XXU*  O.  6661 67uE 
ISJ*  81  XXU*  0.7191232E 
ISJ*  83  XXU*  0.893950  7 E 


ISJ* 
1SJ  = 


85  XXU*  0.906292 IQ 
87  XXU*  0. 1025  32 6E 


rrj*-  P<r  XXU*  <1.117376SE 

ISJ*  91  XXU*  0. 13 7877 JE 
ISJ*  93  XXU*  0. 1667903= 
ISJ*  95  XXU*  9. 2 1 7291 uE 
ISJ*  97  XXU*  9® 1 635760E 
ISJ*  c9  XXU*  0.167624*6 

TTjnorXXU*  0.175729Zc 

UPPER  WAKE  SHAPE  IS 


00  YYU* 
00  YYU* 
00  YYU* 
00  YYU* 
00  YYU* 
OJ  YYU* 
00  YYU* 
00  YYU* 

00  YYU* 
YYU* 

or  YYU* 

01  YYU* 
01  YYU* 

01  YYU* 

02  YYU*' 
02  YYU=- 
02“ YYU* : 


E 0 


0.1919T71F  00 
0. 1985601E  00 
0.1155576E  00 
0.1230558c  90 
0. I 312386E  00- 
0.1 60 3929=  99 
0. 169901 C=  00 
0.1605=586  90 
0.1726313=  00 
0. 18662975  OC 
0.2025991630- 
0.2222686=  90 
0. 2677632E  00 
0.2852830=  39 
0. 2596590= -03 
■0. 1228  659E-0 1 
•0  • 5 1 6 3 R 8 66-07— 


XW( 

1 )* 

0.1757292F  Oi 

YW  ( 

l)=-0. 5160686 E-01 

XW< 

3)* 

0.1663316=  02 

Yw  ( 

3) =-0. 1622 3 1 3E  00 

XW( 

5)* 

0. 220 17635-03 

YVi  ( 

5 1 * 0.13579655-02 

TGA 


SUR=ACF  SHAPE 

0.33M555E-00 

0. 3373662E  00 
9.  3606  661 E 00 
3.3667163E  00 

1)  = 0.1120791=  00 

2) *  l.  63615996-07 

31^-0.5369267^02 

61=  0. 1 361596P-37 

51*  0. I 1332  39= -92 
6)*  0. 19161965-07 
71*  0. I 575 106F-02 
1.1769535F  09  CL30*  U.3521069E 

pontn  * -ano655r2F-oo 

s>r DI2I*  0.1927381=  00 
RCDI 3) = 9.8O89233E-01 
PCD ( 6) * 0*62263 )0=-0l 
■PSSrl  n*  0.2509866=  ul 
■PSSri  21*  0.25617636  01 

■PSSOt  Tl  *-0.26966526-171 

■PSSCI  6 ) * 0.2526677=  Ul 
C070*  u. 

> 9. 1 761 501 E 00  CL  30*  u. 


00 


XXL*  0.2115626=  01  VYL*- 
XXL*  0. 2182256E  01  VYl=- 
XXL  * 0.2269022=  01  VYL*- 
XXL*  0.2315527=  01  VYL*- 
'XXL*  0. 2292566F  01  YYL*"- 
XXL*  0.2650130F  01  YYL-- 
XXL*  0. 25176226  01  YVL*- 
XXL  * 0. 2586835E  01  YYL=- 
XXL*  1.2652  36  7F  01  VYL*- 
XXL*  0.2720013E  01  YYL=- 
'XXC*  Oi2787769=-Or-TYt*- 
XXL*  9. 2855633E  01  YYL  = - 
XXL  * 0.2523599=  01  YYL=- 
XXL * 9 • 2991662E  01  YYL  = - 
XXL * 0.3059816=  01  YYL  = - 
XXL*  0.3 1 28060E  01  YYL=- 
-XXri“0T3T96122T  o r yyc^ 


0. 160258 IF  00 
0.1615332=  00 
0.16266696  00 
0.1635968=  00 
CU16W371 3E-7T0- 
0 • 1669656F  00 
0.1653816E  00 
0.1655959=  09 
9.1656O016  00 
0.1653765=  00 
3 rr6T6  902Y6-0IT 
0.1661680F  00 
0. 1620687E  00 
0. 1615977F  00 
0.15962C96  00 
0.15689626  00 
7TZl5l66C2t  90 


93508  72E- 
3523002  E 


01  CM30*-0. 2362669E  00 

00 


//  SE(  TIMf bS, REGION* |42K> PR r*G# PUN* 1 000  •'  ‘,‘  i , 

rft exee-  portg  ■ — : — * — «--••  - 4—* — •■<■■■ v #-! 

//DECK  OD  8rSOUT*Br0Cfl*f/?ECPM«r83,8lK3IZC*I6»0> 

//COMP  00  ‘ * 

C 3-0  S/C  HYDROFOIL  under  free  surface. 

c PROGRAMMED  BY  fli  PURIIVA  . 3/1/1474; 

C THIS  IS  A PROGRAM  IN  WHICH  OLD  DATA  SXSI(I.j)  ANO  CIRCO(I) 

■t ARE  READ-  FROM  -DATA  CAROS. — 7/4/I9T4. 

DIMENSION  CMM(IO), AGC(IO) 

DIMENSION  WXU(401),WYU(«01),CIR(I0),CIRI(10) 

DIMENSION  AG(lO),Pf>  80(10), CIRCI(?0),BETA0 (100) 

DIMENSION  YBE(3),X7(5$»SXSIO(5,?0),CIRCN(20),CCD0(t0)»ACD(t0) 
DIMENSION  3XSH5),  XXX  (513  ),CP(513), INT( 1 0 ) , XCP (5) , YCP(5> 

DIMENSION  C’lOEI  (200^r«I»E*M»>0)  » XXU  < 1 0 1 ) , Y YU  ( 1 0 1 ) ,-XXl  (UU 

DIMENSION  YYL(10l),AXU(10t  ), A YIK  10 t),AXL( 10 t),AYL< 101) 

DIMENSION  FL(200),FQ(200),FC(7)*PIN(50»7)»OIN(50,7) 

COMMON  FLAPANjCLD, Cl RCD(20),HMH,ALFAZ, SIGMA, SBEIA,XXM,ICPI 
COMMON  IPUl, XA,XB,XC, TANG, FP,YC,YR,UBTG3,XLBIGR, BIGS, 3MAL3,3SS 
COMMON  XSN(5),CCCl,CLE,ERC#YYY'f  XM, ITER A, 3 XSI 0(5), 3X3100(5), VXS(3) 

COMMON  R3I7,LP,SARC(5l34,3AllCn(-513),LPM,0E,I3P,A3Pl,IIII 

COMMON  BETAN(513),HETAM(5l3),IJ,lPK,XII<200),XJJ(200),XOX 
COMMON  A0(7),MPM,MPK , R7ER0D (20) ,CIRCOI (20),NI3P,NNI8P 
COMMON  IFLAG,SXX44 
C PIERO  IS  IN  EON.  14-3, 

C X34LA  IS  A LIMITING  VALUE  FOR  X(«)-X(3)  TO  BE  COMBINED. 

C-MDIV  -IS-USEO-FOR  CALCULAMON-OP-CAVITY  SHAPE -AS  A-  LAST  INCREMENT, 

C NDIVL  0EHIITE3  NDIVI3I0N  OF  LARGER  INCREMENTS  FOR  CAVITY  SHAPE  CALC. 

C NDIVL  MUST  BE  ALWAYS  AN  EVEN  NUMBER. 

C NCAVIT  IS  A NUMBER  OP  INCREMENT  ON  CAVITY*  ALWAYS  MUST  BE  AN  EVEN  NUMBER 
READ(5, 1151)  X34LA,N0IV, NDIVL, NCAVIT 

1151  PORMAT(EI4.7,3I14) 

1 P L A G«  G 

READ (5, 1 1 00 ) NISP,NNI3P 

1100  F0RMAT<2I«)  

C HISP«0  POR  THE  CHANGF  OF  THE  LAST  C0N1R0L  POINT  OVER  THE  SPAN, 
c MNI3P«0  FOR  THE  INCRFMFNT  OP  THF  POSH  ION  OVER  THE  3PAN. 

C 3IGINC  POR  THE  INCREMENT  OP  SIGMA; 

IP(NISP.GF,1>-G0-I0-H4D 

WRITE (6, till,  NNISI* 

till  FORMAT(20X,49M*****LAST  CONTROL  POINT  ON  SPAN  I34****PAI/8*PAI«rl2 
X , 3M/64) 

1110  CONTINUE  

C MCON  IS  USED  POR  CONTROLLING  THE  NUMBER  OF  ITER,  AT  EACH  POSITION, 

REAO(5r2»|-»-NTTCRrM*-MUErN/U-tM<MK,KSTOP,MCON 

201  FORMAT((,IB) 

READ(5,202)  ALFATS,  MHH8,0E, SIGNS, SIGINC  

20 2 PORHAT(SE 14.7) 

C EP3W  IS  USED  FOR  CONTROLLING  OF  DIFFERENCE  BETWEEN  V(IN)  t V(0Ul3, 

C —.USUALLY  FPSQw ASP  1**2.— — 

READ(5r2R3»-AnPr 

203  PORMAT(5E|4:7) 

C PGAP  AMD  FFGAP  ARE  THE  INCREMENTS  OP  INTEGRATION  FOR  FREe  SURFACE, 

C PLI  IS  USED  FOP  THF  LIMIT  BETWEEN  FGAP  t FFGAP, 

C MONK  IS  THE  INDEX  FOR  STOPPING  THE  CALCULATION  OF  THE  FREE  SURFACE 
READ(5*  334)  XLIMIT,FLI 

— 339-P0RMAT(2E 14.7) 

RF AO (5, 344)  NPSUR1, NFSUR2, MONK, NG01,NG02, NCOS 

349  FORMAT (<>IS) 

C NG01»0  MEANS  THE  CALCULATION  OF  FREE  SURFACE  SHAPE, 

C NG02«0  MEANS  PRINTING  OUT  OF  CP. 

C NG03a«  MEANS  PRINTING  OUT  OF  RriP,XMA,XOOtETC, 


FFGAP- 


- KlPM=5  . , 

MPMlaNPN-1 

C MPMI  NO, OP  DISCRETE  POINTS  OVER  SPAN, 

C NITERi  no,  op  calculations  fot  different  stts  of  paraneiers 

C MSTOPI  NO,  IIP  HERS.  WHICH  STOPS  THE  LOOP, 

C HSOPa|  FOR  FLA T»PLATB- CASE, 

c MAX  IT  t NOS.  OF  ITERNS.  WITHIN  TiXFNEH, 

C HHKI  CONTROL  INDEX  FOR  PARAHETERSO 
C NHK>|  FOR  ALFA | 

C NMK»2  FOR  XOLSfMAX.  LENGTH  FO  CHOROi, 

C NHKsJ  FOR  HHHS(SURMER6ENCE), 


NHKsR-FOR  SIGM 


C DE«l,E-4  usually  ih  OXFNFW 
C ASPl  ASPECT  ratio. 

PAI«5,  141592653 
OO  999  IJKL*t. NITER 
ALFAZDsALFAZS 


SIOMAbSIGMs 


HHHaHHHS 

IPTNHK.SC.’tl  RO  TO  240 
IPTNHK.EQ.3*  GO  TO  242 
SIGMAaSlGMStStfilNCwFlOAraJKL-O 
OO  TO  243 

-242  HNHbHMH940.25*FL0AICIJKL-|4 

GO  TO  243 

240  ALP  AZDa  ALP  A ZS-2,  AFLOAT  CXJKL«»I  ) 
243  CONTINUE 
FLAPANaO. 

XMaXXM 

WRITE(6, 998*- ALFAZIb  HHM 

998  FORHATtlX,20H 


ALFA  ZERO»rEl4.7,tOH 


C A3PI  IS  EPSIRONi  l/AR, 

ASPI»l,/ASP  

C RIGS  IS  CHORD  LENGTH  IN  fX«Y),  ALWAYS  2 FOR  THE  RECTANGULAR  PLANPORM 

C-  RIG— I S-  CHlTRO-l  ENG-T  H-IN-fSMALL— XtSHALL  -Y  >-#  A8ED-  ON-SR AN»2, 

C CIRCULATION  TS  GIVEN  BY  THE  DATA  OP  PREVIOUS  CALUCULATION, 

- OO  198  JKOs  1 1 4 

REACTS, 197T  CIRCNTJKOi 
198  WRlTETb.lOSO)  JKU,C!RCN(JKO) 

1050  FORM AT (3X. 14HINITIAL  CIRCNC , 12, 2H)*,El4,7) 

1 97-FORM  AT(E  1 4 .7> 

CALL  COEF (C1RCN, AG) 

CALL  GAHI(CIRCI,AG)  

OO  142  IE« I , KPM 1 

CIRCOITIEtaCIRCITIEl ••••  

142  CIRCD(IE)*CIRCN(IE) 

8TOLL»2»E-4 

8TOLSaR,E«4 

ERC»l.E-2  

CLEal.e-4 

- c CAVIT.  NO.aSIGMA,  AND  P8I7, - 

ALFAZaALFAZD*PAI/lDO, 

WHITE (b, Si  t-V-flIGMA 

511  FORMAT ( lOXjllHCAVIT.  NO  a,E14,7) 

CCC|4AL0G(i.49IGNA)/C2.4PAtV 

C THIS  IS  HTE  8WEPT-SACK  CASE  USED  FOR  THE  3*0  EXPERIMENT, 

BIGO»4,/7 i*4i^A8P 
P3IZ*BIGO/A8PI 

HHM»HHH»8IG0 

OO  615  INPal,MPMl 
REA0(5,61bi  PSSDTINPj 
615  WHITET6,1040)  INP.PSSDTINP) 

,1040  FORMAT (3X,5HPni7<fI2,2Hla,E14, 7 J • . r.( 

61*  format te i a . 7) 


00  000  XIXtat.KSroP 
00  688  ISPal ,HPM1 

.'•••*  Z7HaPAI*(5.-FL0AT(I8P))/8,  ' ' i * 

XF(13P,EQ.4)  0,0  TO  1130  ' 

" 00  TO  list  ....... 

1130  IP<Nl8P.eQ.0>  ZTHbPAI/8.»PAI*FL0AT(NNISP)/64, 

0 Ij-Sl-  CONTINUE 

BIGaB  Jflrt 

•  CZTHaCflSf  ZTH) 

IF(CZ?H.GT.0.5)  BIObBIG-  BIG*(e7TH-0,5) 

' BIG3aftIG*A3P 

c DOWN  IS  OllHN  MASH  VELOCITY  IN  TERH3  OF  ANGLE, 

< OOWNbq . 

DO  148  I0nalr7 

’ 146  DOWNBOIlNN*PLOAT(IOO)»AG(IDnW3IN(IDOa7TM)  

DllWNao.5*DOWN/3IN(ZTH) 

• OOMOamlMN*inO./PAI 
MHITE (6,228)  DllMI> 

•  226— FORMAT (5X, 26HD0WN— WASH— ANGLE — IN-  DEGREE", E 14,7) 

SBETAa-ALFAZ+DOWN 
“ 00  688  ZSEal , NCQN 

WRITE(6,505)  tSP, ISE 

» 505  PORMAT(30X,?5H- POSITION  OF  SPAN**-— , 12,  1 X, IJHNO,  OF  ITER. a, (2) 

WRITE (0, 504 ) BIGS, ISP, ASP 

■» 504  FORNA7(10X,5HHIC3a,Cl4.-VlX,*Ht3PANa,T2,lX»3MAIU,£lftrn 

STOI.ai.E-5 

■>  LPH«7  1 

LPK"20 

» lphi»lfh-i 

LPM2»LFM-E 


LPN4aLPM-4 

c ICP1  13  llseo  FOR  CONTROLLING  PROGRAHi  0 FOR  ITER,  1 FOR  THE  REST, 
C FIND  XSIB.X3IC,X3I0,X8IJ,  USING  • NEM10N'. 

C 3XSI ( 1 laXSIB  

C 3XSl(2laX3lC 

-CSXSl(4TaXSlO 

C 3XSI C41*X3IJ 

ITERAal  

IP(IIII.GE;2)  0.0  TO  744 
IF ( I3E .GE  J 2)  0.0  TO  744 

C DATA  FOR  3X31(11 l NEED  FOUR  CARDS  WHICH  COVER  ALL  SPAN, 

flEAD(5,<,20>-  RX8I(  11^X31(24 r»XSI( 3), 3X31  (41 

620  F0RMAT(4E14.7) 

GO  TO  811 

744  OO  745  LAst,4 

745  SX3I(LA)B3x3I0(LAfI8P} 

8M  P3IZ»PS3D(T8P) 

WRITE  (6,5481-PSU 

540  FORMAT(5X,5HP3lZa,E14.7) 

160  ICPI«0  - 

WRITE (6, 1 02)  ITERA 
102  FORMAT(10X,14HITERATION  NO, a, 12) 

IF ( ITERA .GE.2)  STOLaSTOLa 

IF  (-ITERA  , EO.  M370H1-  MOLaSTOU 

CALL  nXFNrw(3X3l,3rOL»MAXIT>ITN,X34lAl 

630  CONTINUE  - 

231  DO  754  LBai,4 

754  3X310 (LB, ISP) asXSKLR)  

235  DO  537  in|al,4 

X3N(  Illl ) 33X3 1 (-liJt) 

537  WRITE{6,536)  101,3X31(101) 

536  FURMAT(t0X,5H3X3I(,II,?Hla,E14.7) 

SX4*3X3J (4) 

SXX44aSX4 

IFdPlAG.EQ.l)  3X3I(4)bSX3I(3) 


;1  _ ...  . ....  ■ ...  . . i.;  .,kL,.^. 


Iff  !l'L*e.P0.l)*X8NC«)»X8H«»>-  -•  - - 

ICPI*1 

c calculation  or  camma—- circdcispJ 

VA*(SX3I(1 Wl.)*0,5 
VB*tl.-SX3J(l>>*0,5 

ux*sxsiC2Wvo 

UXYallX/VA 

UXY2*UXV**2 

NA*3QRT(UXY2-1.) 

Q2aCCCl*f-VA*WAtvB*AL00(AB3<UXY»WA)J) 

R2«»Vra*3BETA 

WHayB/VA 

32*V A»  30RT(t.-W»***>-VB*(«iB*PAI-AR3TN<HBi» 

UYa3X3lC3>tV0 
IIXXallY/VA 
IIXX2*DXX**2 
MRaSORT  (IIXX2-1>) 

T2*.CCCt«»(  VAaX8»VB*AtnC(UXXfWBjJ 


R2CROa02»R2t32*T2 


CIRCN(l8Pi*R7ERr)*Af5PI*Pai7  /flXSI(A) 

WRITE (6, 515)  I3P,CIRCN<I8P) 

515  F0RMAT(1X,1  tH— — CIRCN(,I2,2HW,etA.7) 

C3PACEa(1  ;«3X31(i))/PinATaPK) 

HC3PAC*0 .5*C SPACE 

FaPACE*C8PAeE/FL0AT<LPH-L-PK-) 

HFSPAC=0.5»F3PACE 

XBET«-l.»C8PACE*FL0AT<LPK-l) 

ICPI»1 

C ICPIaO  FOR  FINDING  8X3T(l),  T.E.,  3X3T  ( I)»YX3(  I)  I ICPl*l  FOR  THE  REST 
C CALCULATION  OF  PRE3SIIRF  DISTRIBUTION  'CP', 

00-25  LP*1,IPN 

C FIND  CPtXSlPl  NEXT; 

SPACEaCSPACe 
IF(LP.f.r.lPK)  GO  TO  91 
XSINa-l . fSPACEoFLOAT (LP-1 J 
GO  TO  92 

91-8P ACE »F SPACE 

X3IN»XBET+  3PACE*FL0AT(LP-UPK) 

92  IMaLP 

IF (LF.EG. 1 ) GO  TO  52 
IF(LP.EO.tPM)  GO  TO  52 
CALL  CCCPPP<X3INf02fIMi 


52  02*1 •♦SIGMA 

53  CP(LP1»1. -02/(1. ♦8IGMAi 

25  CONTINUE 

00  100  LP«I#LPM  — 

FINO  XXXIXSIP)  FIRST. 

XXX(LP)»8ARC(LM 

BETAN(LF»»3BETA 

BETAO(LPla.ALFAZ  

IF (NGM2.GE.  1 ) GO  10  1080 

WRTTE(6, 101)  LP,3ARC(LP>»XXX(LPJ#CP(LP)*BETAN(LP) 

101  P0RMAT(IX»2HI*»I3»lX,5M3ARC*,E14.7,lX,4HXXX»,E14,T»lX»3HCPa»E14,7 

XlX,6HBETANa,E*4,7> 

1080  CONTINUE 
100  CONTINUE 

FINO  LIFT, DRAG  AND  MOMENT  COEFFICIENTS. 

XII<n*-i;/8QRT(l.  + 8IGMA> 

XII(LPM)*SX3I(1 >/SQRT<l. ♦SIGMA) 

POWaSXSI (4  > »PAI 

FSIlaPSII 
OO  105  ITKel,LPM 

IF(ITK.GT.LFK)  GO  TO  106,  , , 

XFSa-l.ACSFACE*FLOATm«*l» ,,  ,,  • ,,  * . , , — , , 

00  TO  105 


tub  xi*amxeet>Ps<»Ace«Pi.o/iT(iTK-iPK) 

105  XJJ<ITK)BPsII/(P0H«(XPS-8XSI(4))i  . . 

Dll  107 

COOET  (IfL)aC03(tlPTA0(IlLM 
8I0ET(ITLi»SIN(8ErAO(ITLU 
roAacp(xTo*xiinrL)*xjj(in) 

PLCTL)»-FBA«CObFT«Tt) 

107  P0(tTL)>F0A*9inETaTL^ 

SPACEbCSPACE 

eLrETa0.5*CSPACf*FL(2)*0.5*F«PACr*FL(LPMl) 

CORAS»0.5*C9PACE*PDCE>*0,5AfnPACE*FD(tPMn 

DO  til  IIIAs2,lPH  3,2 

IP(IUA.GE.LPK»-*PACE»F»PA«! 

CLIPT«CLrPT  + 3PACE*(FL(IUA)+<l.*Pl(IUAflHFL(IUAA2i)/3 
111  C0RAGaCDRARt9PAC£* (FR ( IUA W4.*P0(IUAtl)tF0(IUA+2))/3 
CH»Oi 
IZl»LPH-| 

DO  115  Linal'tZt 

XCEa<XXXfLPM-L10»*XXXC4W»*l3n+liV*0,$ 

XINCaXXXrLPM-LTO»D-XXXaPH-tIO) 

115  CHaCM»0.5*fCPaPM-UnJ+CP<CPM-lI0M*)*XCe*XINC 
CHN(Z9P)aCH 
XLOOaCUFT/CDRAG 
CUIFTaCLtFT/0103 


CORAGaCDRACmr.3 


WRITEC6,U7)  CLIPT,CDRAG,CM,XLOD 

117  PORMAT(10Xt6HCLIPTs,E14,7,lXr6HCDRAGafE14.7f 1 X, RHCMOMENTa.E 14.7, IX 
X,4HL/D»,E14.7> 

CLDaCLIFT*(l.»91GMA) 

CDDaC0RAG*(i;t3IGMAi 

HHITe(6,293>-CL0,CD0 — 

293  rilRMATf  1X,50H- LIFT  ORAG  C0EFFT9  BA3E0  ON  U1 

XIX»4HCDD*,E14.7> 

CCDD(ISP1sCDD*B1G3 

K3TOP1SK3TOP-1  

IF(IUI.LE.K9T0Pli  GO  tO  765 

XP<I3E,b-X;NCCN>-00  £0  743 

IF(NGm.GE.l)  GO  10  765 

C CAVITY  SHAPE,  

OTal.E-5 

XLASTa(3X3I<3»-8X3mi)/N0IV 
GA1JL» (3X91 (31-8X31 <l>-XLA3T)/N0IVt 

GAUS»XLA8f/(NCAVIT-N0IVL4H 

GAL«(-l.»3X3I(2)l/NCAVIT 

NDIallDIVL+2  - - 

NNIaNCAVlTM 

XU»8X3l(li  - 

DO  120  IOHal/NNI 

GAllaGAUL- 

IPn9H.GE.NDI)  GAl|*GAU3 

IF ( I3H.E0.21  XU*3X9X (1 ) 

XI  la  X 1 1 «GAU 

XU»-1 ,-GAl AFLOAT ( I3H-1 i 

IF( I3H.CQ . 1 1 XII«Xt|9XII*DT 

IFrtSW.EO,  1 1-Xl.aXt-ARIMXt-VAOT — 

IF ( X3H.EQ.NNT ) XL»3X8I (2) ♦AR3(9X3I (2) J*OT 
ICPNTaS 

CALL  XINTCG(XUrAXII(I3H)rAYU(I3M)nC0NT) 

ICONTaO 

120  CALL  XINTfO(XL,AXL(I9M),AVUI3M),ICONT) 

AFal./SOOTM  .»9XAMA4 

c CAVITY  3HAPP3  MERE  ARC  RELATIVE  POSITIONS,  THEN  ADO 
C PUR  ABSOLUTE  POSITIONS, 

I0UL*2 
xxut)«o; 

YYL (11*0. 


«4HCL0a,ei4.T, 


XXtt<t»»«i 

YYU{1>«0. 

MhIMallC  AYIT-l 

DO  121  131*1 » NNM»2 

GAMsGAML 

IFtlSI.GE.NDJYU  RMI.CAU8  *AXUCI.SIM>«.AXHaSI*2»*Af /3r- 

XXUCI8I*2>*XXU  »I  ^•**^”isiMWAYU(ISU2»**r/5. 

YYU<I3X*2)»YYJI  131  ISXM>»m(13It2>)*^/S. 

XXUI3I*2>*»XL  191  L i4iUAyL(l9U;,n<Af/5f 

12 1 VYLC  ISI*2 )»YYL(  HI  >-GAl  MATH  4 

WRITPC6/ 12fc)  SHAPE. 2«X.19HL0*ER  CAVITY  SHAPE) 

t26  POnHAT(22X,l«rtllpPEn  CAVITY  SHAPE. 

,23  FaRMAT<9X,/,HT3J*,l3,lX.«HXXW*,El 
XT,lX,flHYYl«.Et4,7) 

X03AaHP8lin, 

03'l«SXSIffl)«-3X3IC3> 

*P  ( JPfcAG.'EO;  1 ) 030*#X3t<*fr 

D|WGsn3/|/PlflAT(ID34) 

TPYD34  LE  XL|HIT)I03««NP9UD2 
tPIDjS.LE.lluiHtTlOING.OIA/rLnATltDSO 

J83«»I03S-I 

OO  130  J3W*l,JOJ4  . 

)H,a3X3I(3)♦i)I,Tfi*F,-0AT(J^,,’ 

,30  CALL^XlMTEGCXtl#  HXlJl  J9H) |WYU(43H)» ICONT) 

xx«cn*xxuaoi> 

YYUaiaVYUtlOO 

101*1034*3  

oa-i3t-tsj*iFtM*^-—  ^^xudsjtn+wxuUSJ+ziT^AP/J, 

til  PORHAftll.lAHUPPEH  BHAPt  IS— — i 

xxui»xxiiuo3«-n 
YYU»*Wtl(X034-l>*P3lz 

XU1»3X3I(<«)tOING 

H,0«10 
-NF3»0 


,36  fSlKcixSsHMH  *<«'*«  **"*] * 

CALL  XlNTEr,(XUl,FXIII,FVU»#ICnNT) 

6CAP»PP6AP 
,35  MP8»HP3*1  . 

_|FlHF8.Ge,MftNK-)-fiO-4R-W 


jPiNPs.eo.t)  on  tom. 

IPCXXU2.LE.FLI)  GGAP»P6AP 

AOS  XII2*Xtll*8GAP  trflNT) 

CALL  XlMTPGf XU2»PXtl2|PVU2# ICON” , 

”5f  ;“5!  «f?5!U.:pyui.ptui>.APAP»A»«. 

XUI»Xll3 
FXU1«FX»3 
PYII1-PYII3 
XXU»*XX+«- 


YYlll»VV>12  . 

IP(HP3.GT>tS) 
iFCxxti2.Gr;io.)  n°  )n 
IF  (XXII2  .GF  .2,  > so  TO  333 
IP <tlP3 .HE J 0 ) GO  Tf).  13^  — 


36 


533  KRlTfifft.lSBl  XXU2,YY U2,NF.3 

1 58  FORMAT(SX,3MXFa,El4.7,5X,3HYFa;E14,7,|X,4HNF3*,I0» 

139  IF(XX||2.GE.«2.)  GU  TO  133  ' 

137  CONTINUB 

, ' 765  WRITE (7, 768  > 8X31 ( I i , 9X31 (2i , 3X31 <3i , 3X4 

768  FORMAT (4E14«7) 

e FIND  THE  NEW  PSII  PnOM-EQN,19-l-4r 

C PIR3T * find  Ron  * THIO, 

‘ C IF  PLAN  form  13  CHANGED,  MUST  CHANGE  THE  ABOVE  EON,  - - 

C XOO  t TOO  ARE  THE  COORDINATES  IN  EONS  29-1  & 2, 

C ROllI  i XOOL , VOtll  ARE  USED  A3  STARTING  VALUES  FOR  INTECRN8.  IN  23-1  ( 2, 
FIHC»|,E-| 

. fin.finc 

_ 220  NV»l 

0 303  TOF-PAI-FIN 

IF(TOP.LE.t.E-l)  GO  TO  302 

* C TOP  DENOTES  THITA  0 PRIMF  IN  P27-2, 

C POP  OEHAOTE8  R-fl-PRIME  IN  26-1  1 27-2, 

, Rnp-(pAi-TnPt*3xsi(4V/8iNtrnP-) — 

WRITE (6, 222)  ROP 

- 222  FORMAT(lX,0HRnP*,E14.7> 

R0P2*RI1P**2 

" 742a3X3I(fl)**2 

CPTrCil3(PAI-TOP) 

» DR0»2.»3X3I(4)*R0R*ei*T 

C POO——  P 27-2. 

■>  C TOO— — P27-2. 

R00a3QRT (ROP2fY42-DROi 

i.  IITaSXSI  (4  )-ROP*CP  T 
T00aARCO3(UT/R0l)i 

„ — c-FINO-IH-t-  RE-IN-EQN8,  25-l-t-  2- — - — — 

CALL  FCNt  (Rflfl,  TOO,  XOOK,  YOOK  , NY) 

- XOllaXnOK 
YOOaYOOK 

» - VOOMa-CIRCD(iaP)/(PAT*Xnoi-OONN 

C MCONai  POR  vCtriCTTY  'V, 

C-MCUNao-  FOR-  aUBMERGENCe— IA4WOI, 

MCONai 

« CALL  OF3IM4(X(IO, VELO,MCON) 

XMAaAB3( VELO-YOOM) 

* IF(NG03.GE.l>  GO  10  1090 

WHITE (6, 304)  XMA,EI*80,XOO,Ynn,VEL0 

" 309  PORHAtIiX,  4HXMAa,E  14^7  rl  Xr6HEP*EP“,E  14.7  rIX^4HXOO*,E-I4,T<  1 X,4MY00> 

XE14.7,SHVELris,Pl4.7i 

a 1090  CONTINUE  

IF(XMa.LE.EPSQ!  GO  TO  302 

n FINaFIN-»FINC  --  - 

NYaNYfl 

COW-  303 

302  CONTINUE 

* MCON«0  

CALL  OF3TMorXflO,YOOn,MCnN^ 

- c HHM  13  FOR  SIIBMt. RGFNCE , - .. - 

206  BOTOH»A3Pl - ( YOO+YOOO) /P3I7 

j.  PSIZaHHH/BOTOM 

P3SD(I3P)*l»3lZa(l,-XXM)tP39D<I3P)#XXM 

» DO  266  NIXnaj,7  

266  WRITE (6, 267  ) NIKO,  AG(NIKClJ 

267  FORMAT(lX,3HAG(,l2,2Hia,f  14.7^ 

C TOTAL  GAMMA  l»RriM  EON,  30-1, 

TCA*,9*AG(  1 ) aPAI 

Dll  272  NI  f AaJ,  7 
XNTl»FLOAT(NlTA-lI 


37  - 


i «ui  •*  1-  » :l> 


r 


XNT2>PI.0At(NI1A*t) 


0 


272  TCAsTUA  + 0.5*AR(NirAW(SIN(XNTt*PAI)/XNT!-3rN(XNr2*P*n/XNT27 
CLAsTSA*ASP/2. 
whiter, 2771  tra,cla 

277  PORMAT(lX,AMTOA«,Et«,7,  1X,$NCL1D«,  E 1 <1  ,7  ) 

ASA  CONTINUE 

C DETERMINE  THE  COEFFICIENTS  OF-FOURieR  EXPANSION  for  circulation,- 
00  376  JllHN*l>4 

176  CIRCO(jnMN)sCIRCN(JnMNl*(!.-XXM) *C t RCD ( JOHN ) AXXM 
DO  327  INK*1,'I 
CIR(IUK)cciRCnf INK) 

127  CIRI(INK)«ciRCOI(INK) 

call  cncrcriRfAC) 

CALI.  GAMItClRIfAO) 

oo  sai  linbi,/i 

WRITE(7,616l  CIRCDaiN) 

1«1  WRITE  (6,5*12)  LIN,CIRCO(LTN> 

5A2  FORMAT (3X,6HCinClH, 1 1 , 2H)  = , P 1 « .7 ) 

“ 00  5*>1  LtN=tf/l 

WRITE (6,61 7)  LIN.PnsoaiN) 

15 1 HRlTf(7,6i6)  P33IHLIN) 

617  F0RMAT(2X,8H— P3SD(,  12,  ?H)  »,  E 1 *i , 7) 

IP(IIII.LT.KSTOP)  f.n  TO  999 

c NEED  A CHANGE  FOR  RIFF,  AR,  A ****** t ************** t *****  * 

C CALCULATE-  3-0  ORAR,-  

CALL  COEF (CCOO , ADD) 

CALL  COEF(CMH,AGC) 

TMs.S*AGC(!)*PAI 
TO«.S*AGDl|)*PAI 
00  1120  N0T»2,7 

XNOl=FLOAT(NOT-iV— 

XN02sPLOAT(N0T*l) 


J 


* ’ 


P01*SINfXN01*PAI)/XNDl-3IN(XN02*PAl)/XND2 
TMsTM  + .5*Ar,CCN0T)*Pni 
1120  TD«T0t.5*AG0(N0T)*P01 

td*td/o , 

TM«TH/R,- 


WRITE(6,!l2lHR,TM 

1121  FORMAT(20X,5HCI)3O*>Elfl,7r3X,9HCH30»,ElA,7) 

C total  Gamma  from  eon.  lo-t. 

TGAa.5»AG(t )*PAI 
DO  279  NIIA=2,7 

XHT1«FL0AT(NITA«4-1 

XHT2»FL0AT(NITAt|) 


279  TGAaTGA+o.5»AG(NITA)*(3IN(XNTl«PAI)/XNTl-SIN(XNr2*PAI)/XNT2) 
CLAsT GAaASP/2. 

WRITE (6, 277 ) TGA,C1.A  

999  CONTI HUE 

STOP 

END 


- 


SUBROUTINE  0P3IM1 ( ANSA,N12) 

DIMENSION  XST (3) 

COMMON  FLAPAN, CLO,CIRC0(20),HHH,ALFA7,SIGMA,SBE rA,XXM,ICPI 
-COMMON  IDOL  ,XAfXB, XC,  1ANC,EP,YC, YR r UBT G.1 , XL3 1 G3, 8 1 G3 , 3MALSf-»S»- 


COMMON  X3N(51,CCCl,CLf.ERC,YVY,XM,nERA,SXSin(5),sXSI0tU5),YX3(5) 
COMMON  P3I7,LP,3ARCC5l3),nARCO(513),LPM*DE, ISP, ASPI,  I HI 
COMMON  BErAN(511),BETAM(5mf  I J , L PK  , X 1 1 C200 > , X J J (200 1 , XOX 
COMMON  AO(7)»NRM,MPK,R78R0D(20),CIRC0t(20),Nl8B,NNl3P  , , 

COMMON  IFl AG , 3XX44 


■ 4 -f  , 

*•  i* 


4o  to  9 

DU  (0  I0«t.4 

10  X3T(IQ)aXSN(I0l 
GO  Til  12 

9 00  11  IH*tf4 

11  X3T(lH)=YX3CIHi 

1?  CSPACEad.PXSTOll/PinAmPtO— 
F3PACE=C3PACC/FLOAT<LPM-tPK) 
LPK2»LPK-2 

xneis-i . ♦esPAce*PioAmi’K-i  * 

ANSA«0. 

UPM331.PM-3 

X3I l »-l .9C3PACP 

BE  1 sBETAN (2l 

APlc(l.tX3n)A(XSr<l)-X3Ili 
Al’lSsSOdT  (API  1 
P3»BE l/APt 3 
Vlia=XST(4l 

iFiNia.EoitV  -mv-ra^ 

IFO.l2.GE.3l  YUSaXOX 
F3*P3/(X31|-YU3) 

5 Oil  1 I«2 , l PM3 » 2 
Fl»P3 

SPACEaCSPACf 

IFd.GE.LPKT  GO  TO  30 

X3I2*-l.»SPACE*PlOAT(ll 

X3I3=X3T2»3PACF 

Gil  TO  31 

30  3PACE»PSPACe 

X3r2»XBET  + 3PACF*PU0AT(t-l,PKM) 
--  X3I3s-X8I2Y3PACE 

31  BE2»BETAN(IM) 

BE3=BETAN(T»2l 

AP?=<l.*XSI2l*(XST(ll-X3T2l 

APT«(1.*XSI31a(XST(11»X8I3) 

A|>2S  = 30RT(AP2) 


-APT 


30RKAP3) 


F2=BE2/AP23 

F3sBE3/AP33 
IP(N12.E0;ii  GO  TO  6 
P2»P2/(X3I2-YU3) 

PTaP3/(X3l3»VU3) 

-fc-F8UMa(P  1 14  . *F2AF3T*8PAC8y-3« 

ANSAaANSAtFSUM 
l CONTINUE 

3QlrSQRr(l.»X3T(tll 

AHTiaBETANf 1 1 *2 . «3QRT (C3PACP ) /SO  I 

ANT2=BETAN(UPN)*2,«30RT(P3PACE)/3Q1 

1 F < N 1 2 . E 0 4 1 V -GO  T 0—7 

ANTl=ANrt/(-l.-YU3l 
ANT2sANr2/fXSTm-YII8i  - ----- 
7 AN3A®AN3A»ANTl*ANT2 

RETURN  

END 


SUBROUTINE  riP3IM2(ANS2) 

DIMENSION  X(3),DA<3l,0B<3l,0C<l>.XTT(3),YYm,THA(3),THB(3> 
DIMENSION  XITB(3l,XITC(3l,Ul(3l«EXU(3),DD(31,0Q(3),M(3),FCN3(31 
DIMENSION  XlTAC31,XST(5),XTrM0513), XITN(513> 

COMMON  FLAPAN,CLP,CIRCn(20l.HMH,ALFA7,3IGMA,3BCrA,XXM,ICPI 

COMMON- -IDlIl  ,XA,XB,XC,  1 ANG  ,EP  , YC  , YR  , UB I GS  , XCB  IGfl  , B IGS«  8MA13  ,333 

COMMON  X3N (S 1 » CCC  1,  d.E, ERC, Y Y Y, XM, II  ERA, SXSIO(S), 3X3100(51, YX3 (5) 
COMMON  PSI7,UP»SARC(313),SARCn(51 J ) ,l.PM,  OE.  ISP,  ASPI  , 1 1 1 1 
COMMON  BElANCjl3),BETAM(5l3),TJ,LPK,XII(200),XJJ(2001,XOX 
COMMON  AO(7),MI»M,MPK,R7ennD(20>,CTRC0I(20),NI3P,NNl3P 
COMMON  IFl AG, 3XX44 


39  - 


on  is  i6>s }* 

13  X3T(r<.)3VXS<T<.i 
PAI*3. 141592653 
C3PACEe(l.»XSrfm/PL0AT(LPK) 
HC3PAC»0.5»C3PACP 
P$P AC Esc SPACE /PLOAT(LPM»lPKj 
— HF3PAC«0.5*F3PACF  


O 


XBET=.i.»C3PACE*Fl.0ATaPK-lJ 

PU*2.«3onr(X3T(m 

PT=X3T(l)-t. 

P3*2.«X3T<|1 

PtJ*PAI*X3T(0l*(l.*X8T(l)l 

PCH3(3)«—  X3TM)/<3QRY(hi«>3ICMA)*XST(<l>*PAI*<XaT('l>-X3T(l>>»- 

FCN3<3>aPCN3(1>*P3IZ 
LPK I*LPM»LPKA 1 
DO  1 IPa|,tPM 

if(ip.eo;m  cn  to  2 

H3PACESHF3PAC 
-3PACE*P3PACf- 


IF(lP.fiT.LPKI)  OO  TO  30 
x(i)BX8rci).sPAce*FLnATap-2) 
X(2)sX(l)»H SPACE 
X(3)«X(1)-3PACE 
GO  TO  31 

-30  H3PACe«HeflPAC 

3PACEaCSPACE 


X(l)»XBET-3PACF*FL0AT(IP-LPKI-l) 
X(2>aX(li-MSPACe 
X(31»X<1^-3PACE 
31  PCN3<l)»FCN3f3> 

NK»3 


IP(IP.EO.IPM)  NK«2 
OO  8 Jo2, NK 

CALL  THlTH2CTHA(l),THB(li,X(l)) 
DA(I)»X(I)»1. 

OB(n»XST(n-X(ll 
-DC  m O30P  T (DA  ( I)  01)0  

xmn>o; 


5 llim»-CCCI*T2.*PAT-THA(li-THB(l) j + XlTTI) 
exii(n«EXP(uinn 
D0(l>oPL*f>C(mxm*PT*P8 
oo(n»(X3T(«ii-xmiAPu 

H(I)»PSlZ*n0TI)/08«> 

8 FCM3m»EXU(n*W(l) 

IF(IP.EO.LPM>  '60  ro  20 
GO  TO  21 

20  FF3«  l./(30RTO>81GHA)«X8r(«)*PAT*(X3T(«l  + l,>) 
FF3*PF3*P3IZ 
-FCMit3ToFF3- 


21  3NMa (FCN3I I >fFCN3{2)*4,AFCN3(31)  AH3PACE/3, 
AM82»AN32*S0M 

IFdO.EQ.Ol  SARC(LPH-IPH)»AN32  . 

60  TO  1 • 

2 3ARC<LPH>«0. 

-AN82»0»- 


t continue 

RETURN 

ENO 


SUBROUTINE  riF3IN«(X0,Yn.HCnNi 
-DIMENSION— V00( 01 rY00J*T0V- 


COMMON  FLAPAN, CLP, Cl RCD(?o1»HHH,ALFA7, SIGMA, 38FTA,XXM,ICPI 
COMMON  I PUL ,XA,Xfl,XC,TANG,FP,VC,VR,U8IG3f  XL  BICS, 0 IGS#  8MAL8,S3S 
COMMON  X3N(S),CCCl,CLE,ERCrYYY,XM,nERA,SXRIn(B),SX3IOn(S),TX3(3> 
COMMON  P3IZ,LP,3ARCC3l3),3ARCn(5l 3)  ,LPM, HE , ISP, ASP! , 1 1 1 1 
COMMON  BEIAN(513),BETAM(?U1,IJ,LPK,XTI(200),XJJ(2001, XflX 


v j V <f  r - v? 


C^MMOtl  ADt71,HCH,M*JK,NZC«00(20)>eifieDI<i*),NlSP,M»»I&f‘ 

: COMMON  IPLAC.3 XX44 

fcl3P  IPENTIPIE9  THE  3PANMIBE  POSITION, 

PAI«J.  141592**3 

YO»0«  

vom«o, 

2"C03(PA|*Pt.nAT(5-I8P)M4V 

IP t I3P.EQ.4 ) fin  TO  30 

00  TO  It  

30  IPtNISP. 6E.il  fiO  TO  II 

2«CO3(PAT/8.»PAI*FL0AT (NNI8PT/64.) 

It  CONTINUE 

ZZaO, 


IPtlSP.eO.n  fiO  TO  4 

ATOa3QRT(XOA*2»ZA«2) 

YTaCIRCOO  TAtXOAATOt 

YIbVT/Za*2 

YXM«VX 

IFCMCn».E9.|}-YI>VJyA*0 

YIM«V| 

OO  TO  5 

4 vis-cinci>(i*/(2,*xoi 
XP(McnN.EO.i)  Yi»cincD(iW(2.*xn*»2) 

YIM-YI 

—5- CONTINUE 

00  1 LAa|,4 
LB»l.Afl 

Z2PaCll3(PAI  AFLOAT  (4-LA  J/«.i 
IF CL*.E0,3T  CO  TO  20 
60  TO  2| 

20  -IP (NI3P.ee  . 1 V-60  TO-2T — 

Z2P«C03rPAI/».+PAI*FLnAT(NNISpJ/64;& 

21  CONTtNUE 
ZlP»(Z2P-ZZW0^5AZZ 
Z2N«-Z2P 
ZlNa-ZlP 

TZ«Z2P 

AiPe3QRrrxn»*2+rz-ziP  S**zS 

A|Ne30PfrXOA*2*(Z-?lN  )*A2> 
A2Pa30PTrxn*»2*r7-Z2PW»2i 
A2Hb3QRT  rX0**2*  t7-Z2N) A*2i 
BlPaXll+AlP 

B|NaXn*AIN- 

i2P«XlHA2P 

B2NaXn*A2N 

IFrMCnN.EQ.lt  BlPaBlP/AlP 
XP(MCON ,E0. I T BlNaBlM/AlN 
IPfMCON.EO. 1 T B2PaB2P/A2P 

IF (MCON.EO* I V-B2NaB2N/A2N 

OAlaCIPCOirLAt 
YOIbOAIaBIP/  IZlP-Z J*a2 
VOM|sGAI«9|N/(7|N-Z)a*2 

IF(LA.E0.41  CO  TO  10 

OAOaClRCD(lB) 

YOM2sCAO*B2N/  «2N  -XWa3 

IPrLQ.C0.I3Pt  CO  TO  2 

Vn?a6A0*B2P/rZ2P-Z)*«? 

60  TO  3 

to  voa«o. 

YOM2ao. 


60— TO- 


2 YO2s-CIPCD(L0)/r2.AXnJ 
IPrMCON.EOllt  YO2=CinC0(L0T/r2.»XO**2) 

3 YOaYOt rvt  + 4. * YOl *Y02T a (72P-7IP »/3, 
YOMaYOMAY YIMA4 . AYOMl aYOM?) A (72P»7 IP)  73 
YI«Y02 


i yrm-  Vo»ti  ’ ■ • — 

nia(V(Hrnn)/<2;*PAr) 
RETURN 
END 


SUBROUTINE  THlTH?(THl,TH?,Xl 

DIMENSION  X3T(3)  

COMMON  rLAPAN,CLn,CIRCn(2Ol,HHH,ALFA?,3IGNA,30erA,XXM,ICPI 
COMMON  IDIJL»XArXBrXC«TANG»EPryC>VR,UBIG3,XLBIG3,BIGSrSMAL 3*333 
COMMON  X3N(5),CCCI,CLe#ERC»trY,XM,ITeRA,SXSI0(S),SXSI00{S),TXS{5V 
COMMON  P3I7,LP#3ARC(5l3),3AUCn(5IS),LPM,0e.I3P,A3PI,  III  I 
COMMON  BETAN(*im,BI?TAM(5t3i,I.M.PK,XtI(200),XjJ(200)#XaX 

COMMON  AOt7),MPM, MPK*R7EROI)(20>, CIRCDI <20 )»NISP,NNI8P 

COMMON  IEl  AG,3XX<I4 

PAI=3. 141542653  

if<icpi.eq;q)  GO  TO  10 

DO  20  13=1,0 
20  X3T(IS)bXSN(1SI 

GOTOSO 

10  DO  00  JT> 1 , 4 
40  X3T(JT>=YX3(JT3 
30  Dt*XST<2)»l. 

02»X3T(2>-X3T(1> 

03»Xtl, 

— 04«X3T<lUx-  

01=01*04 

02=02*03 

A 1 =2 t *80RT (Gl *02) 

Bl=-Cl+G2 

THlsATANCAt/Dt) 

IPHMAT.O;)  THl=P«^T+fl 

H1=XST<3)*1. 

M2«X3T(3)-X3T<1) 

63«HI*D4 

G4=H2*03 

A2=2,*30RT(G3*G4) 

B2»-0i+G4 

TM?=ATAN( A?/02) 

IP(B2.LT;o.)  TH2=PAI*TH2 
RETURN 

ENO  ..  . 


SUBROUTINE- GCCPPP(X8I,Q2r IS) 

DIMENSION  INT< 10) 

COMMON  FLAPAN,CLn,CIRC0(20),MHH,AlFA7,3IGMA,3BErA,XXMf ICPI 
COMMON  l0tlLrXA>XB,XC,TANG,EP,YC,VR,UBlG3,XLBIG3,6IG3,SMAl.9>333 
COMMON  X3N<5),CCCl»CLf,ERC,TTT,XM, UEIJA,  3X810(5  >,  SXSIO0(S)  , YX3<5> 
COMMON  P3IZ,LP.3ARC<Sl3),SARCnT3lI),LPM,0E,I3P,A3PT#IIII 

COMMON-BET  AN<*il3)rBETAMC^mrMrLPK,  XU  (200)  rXJJ  (200  ),XOX 

COMMON  AO(7 ) * MPM, MPK ,1)71:0(10(20 ), CIRC 01(20), NISP,NNI3P 

COMMON  IPL AG*  3XX44  

JI8=0 

PAIbJ. 141542653  

ACPl=30nT((l.»X8I)*(X3N(l)-XSli) 

ACP2=2,*0GRT(X3NU->4 

ACP3»ACP2*ACP1*XSI*<X3NC I Wl . )»2, *X3N( I ) 

ACP4»1.*X3N(1 ) 

ACPS»ACP3/ACP4 

CALL  TM1TM2<TM|,TH2,X3I) 

1 P3=0» 

— 2-WAb-CCCI  »(2,  *PAI-TM4-TMi)>BLj 

EUAaEXP(IIA) 

XII(IB)bEUA*ACP5  

Qb-xSI/XIICIB) 

, G2«G**2  - f . - : 

RETURN 


SUBROUTINE  OXPNEH ( X, ATHL , M, I, X 341  A)  ' 

OOUBLE  PRECISION  A0PlfA0P2,ABP4,00P| 

DIMENSION  F(4),P(50,4),X(4),INT(l0),G(4,4) 

COHMOll  FLAPAN,Cin,CIRC0<20),HNH,ALFA7, SIGMA, 3BCfA,XXM,ICPl 

COMMON  IDOL, XA, XB,XC,1ANC,fPrTC,TR,U»IG9,XUMGS, BIGS, SMALS, SSI 

COMMON  X3N(5>,CCCl,CLk'.eRC,  YVr,XM,  HERA, SXSIOtS ) , 3X3100(5) » YXS<53 
COMMON  P3IZ,LP,SARC(5m,SARCn<SI3),tPM,0E,I3P.ASPI,IIII 
COMMON  8ETAN(S13),BETAM(5|  3 > , I J,IPK , XI I (200 ) , X JJ (200) , XOX 
COMMON  A0(7 ), MPM, MPK, R7BROO{20>, C IRCOt (20), NIBP,NNI8P 
COMMON  IFL AG, 3XX44 

IPLAGaO 

CIRC*CIRCD(13P) 

IZZaO  

0G«l.E*3 

df*i.e-5 

00-67— II  Jaf-yO * 


67  HRTTE(6,66)  IIJ,X(IIJ$ 

66  FORIlAT(iX,2HXf,II,2H)>,et<l.7) 

55  SI 1B2. *DE 

SI2a-l.a2.40G 

IP(X<n.LT,8ll)  X C I ) *3T 1 

IP(X(2),GT.SI2)  X(2)*SU 

112  IFCX(2>.LT.-9.)  X(2)a-4. 

B10aX(n*2.*DE*X(l> 

IP(X<3).IT.3M>  X(3)«3I4 
IF(X<3),CE.5.*>  X(3)*5, 

SI5aX(3)42.«Dr*X(3) 

e-  IF( X (4  ) ,U  T»UE-4  V -X  ( 4 )ax  (-3U-4- X(4  ) a X ( « ) -X  ( 3)  ARE-  USED 
X34AaAB3(X(4)-X(3>) 

IF(X34A.UT.X34UA)  GO  TO  171 
IF(X(4).UT.O.)  GO  TO  1055 
IF(X(4).LT:x34LA)  GO  TO  169 
IFdFLAG.EQ.l)  GO  TO  1065 

IF(X (4V.L-T JSI'it  -Xd'aSTS 

GO  TO  170 

1055  IFUAGal  

X (4)  axonal .2-2 
GO  TO  170 
169  IFUAGal 

CO  TO  170 

1065  IFLAOaO 

X(4)ax<3UX(4> 

GO  TO  170 

1 7 1 X (4  )*X34A  

IFUAGal 

170- COOT  I WOE 

DO  6(1  1 1, I*  | , 4 

• 68  MRITE(6,66)  IU,X(IIJi 

P*I«3. 141592653 

IJa  l 

C FIND  F( 1 ) . 

00-20— IKa  |^4 

20  YX3(IK)aX(IKt 

5 AiarXBdUt,  

A2aVX3(2Wt; 

A3ayX3(3)4|.  

A4*VX3(2)-YX3(1) 

A'ia  yX3(3  V-YXA(  I-) 

8iacORT(YX3(l)» 

*0P2*A2 

AOPIBAI 

ADP43A4 

00P|a.A0Pt/(2.40SGRT(ADP2aA0P4)4ADP2»A0P4) 


Vi~DOPi  •-*  -r: 

DlLflS=ALnC(01 ) 

02«At/<2;*80RT(A?*A?)»A3»A5i 

D2L0GaALnC(D2) 

' *■  THA|»2.aB| 

THA2*yX3(l)-l. 

ATHA2aAB8CTHA2> 

THAa-ATAN(THAt/ATHA2) 

IF (THA2.LT  ^0  , ) TMAa-PAl-THA 

D3aPAI*THA 

ANSaSneTA 

2 FAa.(cCCI*(DILnG+D2LnG>»0J*AN8i 

iFttj.eo.n  F(t»spA 

IF ( I J .CO .2 ) GO  TO  3 
lPdJ.CQ.3)  CO  TO  4 
P(l,2Ta-CCCl/30RT(A2AA4) 

Pd#3la.CCCl/3QRT(A3*A5> 

IJ»2 

VX3 ( ITaXfi j toe 

00  TO  5 

3 P|Pa-PA 

IJ»3 

vxschaxd  i-oc 

GO  TO  5 

4 FlQr-PA r 

PM»na(FlP-FlO)/C2.*DE) 

P(W4)I0, 
c HEX  T ABOUT  F(2) 

IJa4 

00  30  IMe|,4 

30  YX9«ttb*Xd 

8 AAl=YX3(lWl. 

AA2aYXS(2)«l,  , 

AA3=YXS(3Wi; 

IP(IFLAG,EO.I)  YX8(4*aX(3l 
AA4aYX8(4 )♦! . 

AAr,aYX9(2i-YX9(<T» 

AA6=YX3(4)-YX3(3> 

IPdFLAG.EO.t)  AA6aX(4i 
AA7aVX3(4}.YX3U) 

AA8aYX8(2)-YX3(h  . 

AA9aYX8(3t-YX3d) 

BBU l»A A | *A A3 

BBt«AA4*AA7 

BB?sAA2*AAB  - 

BB3>AA3*AAO 

3Q301a3QRT(BOl> 

BQB82aSQRT (B02) 

8QBB3a80R  F (B03 V 

BBLta2.iiS0RM2*90nfltAYXS(4)«(AA2>AAH)fAAn-YX3(n*AA2 

D01LOsALOG(BQUl/BflLI ) 

BBU28AA1*AA6 

BBL2b2. *3QBB3 aSOBB I *YX3(4  J a ( AA3>AA9 J ♦AAO-YXSd ) AAA3 
002LOaALOG(BBu2/BBL2i 

TH0|aj.»80RTfVXBd>-U«O»*l 

TH32=YX3(1)*AA4-AA7 

ATHB2aAB3(THB2)  - 

THBa-ATAN(THnt/ATHB2i 
IF<THB2.LT;o.)  THBa.PAI-IMB 
PATaPAl+THB 

AU8 | asBET A 

7 FBa-(CCC t a(00 lLOt0O2LO) FPATAAN3I ) 

IFdJ.C0.4)  F(2)aFB 
IFdJ.EO.S)  00  TO  « 

v.m.  iFdJ.ro.o)  on  to  i*..,. »... *-•.  . , • 

TFdJ.EO.7)  GO  TO  II 


lPdJ.Ea.81  60  TO  IE 

iPdJ.eoJm)  r.o  to  427 
IFdJ.EO.102)  60  TO  428 
PPl«l./AA5 

PP2»3Q0B1*(AA2*AAB)/8QBB2*AA4*AA7 

PPJ»PP2/BBll 

— P(2,2>sCCC|*(PPl-PI»T) 

PPOa-t ./AAfa 

PP5a3QBBl*(AA3*AA9)/306B3+AA4+AA7 
PP6aPP5/B8l2 
IFdFLAC.eO.t)  60  T(1  425 
P(2f3)aCCCl*  f PP4-PP6) 

— 0070-426 


425  IJ4101 

VX3(J)«X(3W0F«*X<3) 

TX3(4J«YX3(3) 

60  TO  B 

427  F5P«-FB 

1 J»l  02 

VX3(3)«X(3)"0F*XC3) 

VX3(4)«YX8(3) 

OO  TO  « 

428  FSQa-FB 

P(2,3>«(F5P-FS0)/(2,*0F*X(3)) 

— >426  CONTINUE- 

IJaS 

TX3(nax(|)»l>e 

GO  TO  B 

9 F2Fa«FB 

ij»6 

VX8U)«Xd>-0« 

00  TO  B 

10  F2 Oa.FB 

P<2# l)4(P?P»P20)/(2,*0E) 

U«7  

YX3dl«XO) 

iftiflag;fo:t)-6o-to-ib2 

TX3f4)ax(/|)tOP*X(4‘ 

• - GO  TO  B 

It  F3P»-FB 

1 JaB 

VX3(4)ax  (4>«I)F4X  (4) 

60—70-  B- 

12  P30a»PB 

P<2#4lB(F3P.P3Q>/(2,*0F*X{4)>  .... 

GO  TO  1B3 

182  PC2,4)aCCCl/X(4) 

| 83  CONTINUE 

C-NEXf-FCm 

OO  40  IHB1.4 

40  YXS(IM)aXdM)  

IF(IFLAG.EO.l)  YX3(4)a  X(3) 

VAlaYx3(2)-YX3T4)  - - 

YA2=YX3(3)-YX3(4) 

IF(  IFLAG.E0.-1 1 — YA2*"X <4) : 

YYI»YA1/YA2 

••  YY1I.«A1,0G(YY1)  - 

HPacinC4FAT*YX3(4)/(A3PI*P8I7) 

F(3)a-(vxS(2)-YX3(3)fYYH.*YXS<4)-HP) 

p(3,n»o. 

YA3aVX3(4)/VA-l 

VA4aVX3(4)/VA2 

P(3,2)sl.*YA3 

Pn»3)*-i;-YA4 

IFdFLAG.eO.l)  P(3f3)aAL06C(YXn(S)-YXS<2))/YX8C4))*YX3(2)/ 
1 { YX3(3 )-YX3(2) ) 


— - ir«I^LARirO.I>  60-T0  262 

P(3»fl)»Vm-YAl+tAa-HP/VX3<4) 
60  TO  263 

262  P(3'4)s.YXS(3)/m<«) 

263  CONTINUE 
C NeXT  PC4). 

14610 

vxs<n>x(i  woe 
19  CALL  0F3IM2TAN32) 

iP(iJ.eo,9)  cri  to  7i 
IF(IJ.EQ.IO)  60  TO  13 
IF(I4,EQtll)  60  TO  14 

IPtlJ, eo;»2»  60-70-  15 

IF(IJ.E0.13)  60  TO  16 
IF(IJ.E0;t4l  60  TO  17 
IF ( 14 ,E0, 15)  60  TO  1R 
IF ( 14 ,EQ, 16)  60  TO  21 
IF(IJ.E0,17)  60  TO  22 


— 1TAN3P»ANS2 


IJ«|1 

vx3<n«x(n»oe 

CO  TO  19 

14  AN3Q«ANS2 

P<4,  I )«.(AN3P-ANfl9)/<2.*DE) 

VX3«lT6X(t1-  - 

I J»9 

60  TO  19 
75  AM3F*AN82 

F(4)"»(SIC3“AN9F) 

14612 

VX3C2)«XT2)tl>6*AB3CX«)4 

60  TO  19 

15  AHSPP>ANS2 

14613 

TX3(2)«X(2)-OG*ABotXT2>> 

60  TO  19 

16-AN86Q»AN82 ■ 

P(4,2ia-(AN9PP-AN300)/T2.*OCAAB.1fXf2))> 

VXS(2l«X(2) 

14614 

YX3(3>6X(3)*0F«X(3) 

IF(IFlAG.EQ.I)  VXS(4)6VX8(3) 

GO  TlV  19 

17  AMS 1P*AN92 
IJ615 

TX3(316X(3)-0F*X(3) 

IF(lFLAG.eo;n  TX9<4)6»X8(3) 

60  TO  19 

l*-AN8tQ6ANS? 


P(4.3>6.(AN31P-AN310)/(2.*6F»X(3)) 
VXS(3)6X  (3) 

IP(IPLAG.eQ.l)  60  TO  371 

I4«14  

VXS(416X(4>T0P«X(4) 

«tf~T<M9 

21  ANA6AN32 
I4«17 

VX3(4)»X(4)»0F*X(4) 

00  TO  19 

22  ANB6AH82 


TX3(4)6X(4) 

00  TO  372 

371  P(4,4140, 

372  CONTINUE 

CALL  0ETEBM<P#4#DeT8O) 


90  25  IOETal.4 
DO  26  LPOst .4 

oapo,XDETV»PCLPcfioBTi  > 

26  pclpg,ideiv»pclp«) 
call  i>eTer»M<p,fl,neT«* 
rPCXDET.EO.lV  0ELB«0FTE/0FT6O 

— IPCIDET .E0.2 V DELCaOFTE/DBTBO- 
IPCIDET.E0.3V  DELDBDFlE/DFTOn 
TPCTDET.EO^V  Dri JaOFTE/DETBO 
OO  27  LPG*| , 4 

27  PCLPG,IDETVbQCLPG,IDEtI 

25  CONTINUE 

X40»X(4V 

XC1)«X(1)+0ELG 
XC2)»X(2)+DELC 
XC3)aXC3)«DEL0 
XC4)»X C4)+0EU 
DO  60  LMN*t , 4 

60-WKlTE(6*61  V LNNf X (fcNN-J 

61  P0RHATClX,2HXC«lt»2HVs«E14.7) 
ABSOSABSCDELO/XCIVV 
ABaC«ABsC0ELC/XC2)l 
AB3D«AB3(DELD/XC3VV 
AB3J*ABS(DELJ/X(4)V 

KEIQ«0 

XPCABSB.LT. 8T0L)  KEXfl»l 
XPCABsc.cr.srnLV  keio»o 
IPCABSD.CT.STOL)  KEIflsO 
IPCAB3J.Ct.9TUL)  KEIO*0 
IPCKEXO.E0.lV  CO  TO  33 

I»I+l 

WRITE (6.42)  I 

42  P')RHATC20X,14H1TFRATI0N  ND,«,I! 
IPCX.EO.N)  GO  TO  33 
GO  TO  55  

35  IFCI.EO.MV  CO  TO  36 


GO- T 0-36 — 


36  WRITE C6.37) 

37  FORMAT C 1X» 34H0XPNEW  DIO  NOT  CONVERGE  WITHIN  > M » 3 
IFCXC 1 ) ,LT ,31 1 V X(1)=3I1 

XFCXC2V.GI.BI2)  X(2V»SI2 
IFCXC2V.LI.-9.V  X(2V»-4. 

SI4«XC!>*2,*0E*XC1V • 

IFCXC3)  .LT.9I4)  X(3V>9I4 

IPCXC3V.CF.S.)  X(3j»3, 

9I5»X(3)t?',  *DP*X(3V 
XPCXC4V.LT:o.V  GO  TO  1132 
XFCIFLAC.eO.l)  GO  TO  3B0 

IFC  X C*V,L  1.6 15 )X<  6)*6tS 

360  CONTINUE 

1152  XC4)»X40<»i;e-2  

SB  DU  129  ITX«1,4 
129  WHITEC6. 13| V ITX.FCITXV  - 

1 31  FORMA TC IX.2HF (,tl,2H)«,EI4,7V 

DiviJa-iup.H^ 

on  132  IIIObI'A 

132  WRITEC6. 133)  IUP,IUG,PCXUI»,IlloV 

133  P0RMATCtX,2MPC,Ii,tH,,Il,2H)«,E14.7V 
RETURN 


SUBROUTINE  XINTEG(X,XC,V6, ICONT) 


COMMON  PLAPAN,CLI),CIRC0C?0V,HHH,ALFA7,9IGMA,SBErA,XXM,lCPI 
COMMON  IDOL, XA,XB,XC,TANC,EPf TC,VRfUBI63,XLBICB, BIGS, SMALS'SSS 
COMMON  X3N(5V#CCClfCLE.ERC.TW,XMf TTERA, 3X310(3) ,9X3IOOC5),rXSC3) 
COMMON  P9!7,LP/!IARCr513),3AnCnr5t3),LPW,0E,I3P'A8P!,TIII 
f OR MOM  »*rAM<*l3>,BETAM(5l3>,IJ'LRK  >U  I <300), KJJ <*••> .XOX.  ...  ... 


COMMON  An(T)>Mi*M,Nf>K,n?ERnO<20),CIRCOI(2O)'NIBf>tNNiaM 
COMMON  IPIaG,3XX44 
PAIaJ. 141342*13 

Al>XSN(l)+|.  

A2»X3N(21-X 

A3«X3N(2)*t; 

A4*X3N(2)«XSN(1) 

A3«X+1. 

Ab«X-XSN(l) 

B2«X3II(31-X 

B3»X3N(3)M. 

»4*X3MC31-X3Nm 

IF(ICONT.EQ,4)  A2«-A2 

IPUCONT.EO.X)  82a-B2  

XP(ICONT.EO.S)  B?a-B2 
Cl«Al*A2 

C2*30RT(A3*A4*A5AA<>) 

-CjsX*(A3*A/|V 

C4«A4-XSN<1)«A3 

0l«Al*B2 

02a80RT(B3*B4*A5*A<iJ 
D3»X*  <B3  *84 1 
04«n4-X3N(i  WB3 

0taCl/<2;*C2*C3+C«> 

F|«01/(2,*02*D3*D4> 

31MlBCCCl«(ALflG(Rll*ALnG(Fli) 

IF < ICllNT.EO.4 ) »IN!*-3INI 
AAUsX»(X3N(li-l.)*2.*X3N(l> 

AAl«X«A| 

IF<IC»)NFiE0.4->-AAl.-»»AAt 

BC&aAAl»/AAl 
THP«AUSIN(BCO) 
lP<tenNT;EQ.4>  tnp»-thp 

TH2»0.5»PAI*rMP 
Hla-3QRT(A5‘A»>/P- I 

-fP<IC»Wf-,E0.4>-+*la-M4 

XUXaX 

CALC  OF 3IMI(AN3A,1C0N1> 

3IN2«H1»AN3A 

S3Aa3INl*TH2t3lM2 

C033»Cn3(33A) 

3I33«3IN(33AV 


X3N4»X3N(fl  > 

IF< IFL AC ,PO. | ) XSN4»X3N<3>*3XX44 
3POa-P3lZ«X/(X3N4  *PAla<X-X3N4 

XG«C03S»3PO  

YG*3I33*3PO 

HE TURN 

ENO 


SUBROUTINE  PCNl(P,TOP,Xnx,YnK'NY)  , 

OIMEU3ION  X3TI5)  

COMMON  FLAPAN,CLrt,C|RC0(?0l,MHH,AtFA?,3IGMA,3BErA,XXM,ICPI 

COMMON  If>Ul»XAfX»,XCTl*MC»EPr-YCrYnrUBr63,XCBIGSrBIG3r3MAL3r»»» 

COMMON  X3N(b>,CCCl,CLE.EfiC,YYY,XM,nEnA,8XSlD(S),3X3T00(5),YX3(3» 
COMMON  P3II, UP# SARC(3l3),3ARCn(51S),LPM,0e. ISP , A8PT , I I I I 
COMMON  BETAN(513),BETAM(5mf  U , L PK, XI I (200) , X J J (200) , XOX 
COMMON  AnC7i»MPM,MPK,nrEnnO(20),CrRCOr(20)»NI8P»NNl«P 
COMMON  IFl AG  1 3XX44 

PA  I a)  . 1 4 15«2b*»l 

SlaP/FLOAICNY) 

C FIRST  FIHB  XX|  IN  P25.|,  . 

XOOK«0,  , 

- , YO0K«0,  : • - r-H — - .-rr:-.  -v 

XOKaO i 


VOfta#. 

» oo  to  r*t,NY1  , 

! Rb3I»FI.OAT(I>  1 1 - I*  ' ! 

AsR*CO3<T0P)»  X3N(4) 
BaR*3IN(TQP) 

R2*30UT(A**2T8**2> 

T2*ATAN(8/A)  

IF(A;lE.o;>  T2»PAIM2 
XKI*-PSIZ*R/(PAI*  X8N(4)»R2i 
c NEXT  CALCULATE  IM  $ RE, 
AA*R«cn3(?0P)»t. 
BB*R»C03(TOP>-X3N(t) 

cc*r»8IN( rnP> 

RAa30RT(AA**2+CC**2) 
RB«r>0RT(BB»*2*CC*»2i  .... 

R t *SORT(RA*nn) 

THAaATANfCe/AA* 

TMObATAH(CC/OB) 

ifcaa.le  ;0;S-Tha«pa-i^t«a 

IF(BB.LE,0.  ) 1 HBaPAI+THB 

T1»0«5*(THATTHB) 

XRRaO. 

XMHaO. 

DO  6 MIOat.O 

CALI.  PHI  NT  (AAR,AAM,MI0^R,-WIP-) 

IFfMIQ^EQ. t S XRa-CCC  t A AAR 
IFOIIQ.EO.I ) XMa-CCC t *AAM 
IF (HIQ.CQ.2 ) xnaCCCtaAAR 
IF(MIQ.P(j;?)  XHaCCC  1 *A  AM 
IP<MIQ.EQ.3>  XRa-AAR/PAJ 

IF(WIO.eo:3»-XP«-AA»i/PA-l 

IFCHIQ.EQ.4)  XRx-AAR 
XF(HIQ.EQ.A)  XMa-AAH 
XRRaXRR+XR 

a VHH>yHM4VM 

PIMaRl*(XRRASIN<TO»XMM*Cn3(T|} 

PREaRt»(XRP*CnS<TlV-XMM»SIH<TtV 

EXlsEXP(-PIM) 

ARlaPRE+2.»T0P»T2 

000*A3PI*XKl*eXI 

XOKbOOO*C03(AR1) 

YQK*0DD*3IN(AR1) 

XOKa  <XOOK<M<JK  >*0.5*81f*nR 

yOKa(YOim*Y(JK)aO,5«8t*VnK 
XOOKaxOK 
10  YllOKayOK 

RETURN  

CNO 


BUBRniJTINE  RMTNTtSR, 3H, M1Q, R, TOO) 

DIMENSION  RriT(200), X3T(B),XKERl(2(»0>,XKER2r200t,XKERM(200> 
DIMENSION  XKERRT200) 

COMMON  FLAPAN,CL0,CIRC0<20),HHH,ALFA7,8lGMA»S8ErA,XXM,ICPI 
COMMON  IDIIl  ,XA,XB,XC,TANG,PP,  VC,  YR,  UBtGS, XLBIGS, BIGS, 3MALS, 333 

COMMON  X8N(S*rCCClrCLerEIJC»-YXY,XM,mnA,3XSln<M,3X3I0n<5)fJfXJA 

COMMON  P3tZ#UP»SARi:(515),SARCn(5t5),LPM,0E, ISP, ASP! , I I I I 
COMMON  BETANCil J),nETAM(51ii,IJ»LPK,XTI<200>»XJJ(200>,XOX 
COMMON  An<71,MPM,MI’K,R7ERnD(20>,CIRCOI<20),NI8P,NNISP 
COMMON  IFL  AG,3XX44 
PAIa3.Mt5<»2h,53 

UPMJ=LPM-5 

IF ( ICPI , EO. 0 ) so  ro  to 

00  12  10*1,4 
12  X3T(lO)=X3N<Ifl) 

CO  TO  It 
10  00  t 13=1,4 


1 *$TCTS)*YXSCt})  ' - -~ 

u xxj*R«cns<Tfin» 

YY  taR«SIN (TOO) 

SP|s-l.-X3T(2) 

SP2oXsTm-XSTM> 

SP3sX3T(tUi; 

3P'lnX8T(tV  — 

IFfMIQ.EO.3*  GO  TO  2 
IFfHIQ.EO.tl  3INC*BPt/2t« 
IFfHIQ.EQ.2i  3INC»BP2/?i; 

ip (mio.eoI^)  siNCBr.p«/?i- 

IPfMQ.eQ.il  X3TAR=X3T f 2 » 

IFfMIQ.EO.2i  X3TARoX3f  f H+SINC  - 

IF(MIQ.EQ.A)  XSTARboJ 
OO  3 ITP«l,2l 

X YIN»X 3 TA REFLOAT tITP«|}»3INC 
XIN8eXYlN-X3Tf l> 

IF(HIO.EO.A)  XIN3»-XIN3 

PTaSQRTf  tXYINMi>*XIN3i 

Pll»fXYIN  -XXI ) **2+YY  I **2 

XKeni(lTP)sl./(P7*Pll) 

3 XKER2aTPi  = (XYIN-XXl)nXK2Rl(ITP> 
8R«0. 


OO  \ -Itn*tit<»,2 — 

3R»3R+SINC*(XKFR2(TLn1»4.*XKER?(TLnHi*XKER2(ILfl+2'))/3« 
a 3M*3Mt31NC»(XKeRt(ILO)4-4.*XKERt  ( ILO»  I ) ♦XKt'R  1 (H.fH2n/3, 
SM3SM»YY1 

IF (MIQ. E3j  1 1 AI)Qe2,,*flQRTf3INC)/((f-l.-XX|>**2*YYI<**2)(* 
XSQRTf  i,  + XST<n1> 

- - IP<HIQ.EQ.2>  ADD»2f*flQRTf8INCi,M(fXSTt»1-XX»1*«2*YYI<M2l* 
X8QRTfX3Tfi>M.i) 

IPfMIO.EQ.M  AODe2,A8Qnr(3INC)/(((X3T(I)-XXn**24'YYl*A2)* 
XSQRT<t.*XSTtt)1) 

IF(MIO.EO.l)  ATOPs-l.-XXI 
IFfMIO.eQ.2i  ATOP=XST(|)-XXI 

IF(HlQ.£0;»>--AT0P  = X87<»-V*X*t 

3H*SR»AOO*ATOP 
SMbSH+AOPaYYI 
OO  TO  S 

2 C3PACE*  ( l .+  X3T(m/FUflATaPK) 

F3PACEoC SPACE /FLOAT (LPM-EPKi 

LP«1»LPH-| 

XIITe-li^CSPACFAFLOATaPK-li 

OO  (,  KIK*2,LPM1  

SPACPaCSPACE 

IPfKTK.GfJlPKl  3PACE»F3PACE 
XPl=-t.*3PACE*FLOAT(KIK-l) 

IFfKIK.GT.LPKV  XF|aXHT*Fl-OA.T<tHK-LPK)«3PACe 

IFf  fTERA ,EQ. 1 ) GO  TO  20 
XTOPRaSE  f AN <K IK ) A <XP l-XX 1 i 
XTOPMsBET ANf K IK  1* VY I 

GO  TO  2|  

20  XT0PRo3BCTA»TXFI-XX| J 

XTOPMbSBGTAaYYI 

21  XHOTAaSQRTf fl.AXFti*fX3Tf li-XPll) 
X»0T0»(XFl-XXn*«2*YYI**2 
XB0T«X3nTA*xnilT# 

XKERRfKIKi«xTOPR/X80T 

b XKERM(KIK)*XTnPM/XBOT 


sm*o; 

00  7 KIN»2«LPB3f 2 
3PACG«C3PAC2 

IFfKlM.GEAPKi  8PACE»FBPACS  , 

3R«flR*3PACE*(XKrRRfKTM>»4,»XKEnn(KtM»l )»XKERR(KIH»2i)/1 


50  - 


7 SM=SM+SPACE*(XKERM(KIM)+4.  *XKERM(KIM+l)+XKERM(KIM+2)] 

sibxsqrtci.axstci)) 

SAODtaBeTAN(t) 

BA0D3aBCTAN(UPM) 

ifcitera.eq.i)  badoibsbeta 
ifcitera.eq.i)  baodb«sbeta 

— A00i=BAD0|*2,»80IUCC3PACF)/<»TB*((-I.-XXO*«2*VVI*a2)) 

AODa*BAODB*2.*3UnrfF3PACF)/(RIB*C(XST(n-XXn**2»YYt*»?)» 
3R»TR»A0i>l*(.l,-XXn*ADD8*CX3T(n-XXt  ) 

SM  = SM  + AI>D| * Y Y t ♦ ADDBa Y Y 1 
5 RETURN 
END 


8U0ROIITINF  COEFCG,A) 

c COEFFICIENTS  OP  FOURIER  SINE  8IRIE3.  

DIMENSION  GC4),AC7),PC00,7T,0C50,7i 

COMMON  FLAPAN,CLn,CIRCDC20),HHH,ALFAZ, SIGMA, 3BETA, XXH, ICPI 
COMMON  IDIJL ,XA, XB,XC,TANG,FPfYC,YR,U0IGn,XLBIG8,BIG3,SMALS,333 

COMMON -XSNC5 ), CCC 1 ,CLC. ERC, YV Y, XM, I TERA, 3X810(5) ,SX3I00(5),YXS<53- 

COMMON  PSI7,LP»»ARCC513),nArtCn(5l3',LPM,DE,ISP,ASPI, III  I 
COMMON  BETANC5m,BETAM<5l3),U,LPK,XII(200),  XJJC200),XOX 
COMMON  AO(7),MPM,MPK,R?ERn0(20),CIRCDI(20),NI3P,NNl3P 
COMMON  IFL AG, 3XX44 
PAIb3. 141542653 

Da-4  -I«lr7 

DO  1 J»l,7 

1 P(I,J)«3INCPAIaFLQATCIT  AFLOAT C J ) /3 , T 

IF (NISP.GE  , t ) GO  TO  6 

DO  7 K Na 1 , 7 

P( I, KN)aSI N C CPA 1/B, AFLOAT CNNI3P ) APA t/64,) AFLOAT (KN)) 

7 P (7-r  MHVa8IN(-fPAI»PA4-/TT,~FtT‘A-T-(NNT3P  >APA  I /64,)aFL0A  T (KN)-T — — — 

6 CONTINUE 

CALL  DErERM(P,7,Bi 
DO  2 1=1,7 

DO  3 J=1 , 7 

o(j,n»p(j,i) 

IFC.J.GE.sV  GO- TO  4 

PCJ,I)bGC5-JT 
GO  TO  3 

5 PCJ,I)bGCJ-3) 

3 CONTINUE 

CALL  OETERm (P,7»C) 

AT  I )bC/3 

DO  4 K»  1,7 

4 P(K,n=0(K,lV 

2 CONTINUE 
RETURN 
END 


SUBROUTINE  GAMI (C, A) 

DIMENSION  CC10),A(toT«X8CtoT 

COMMON  FLAPAN,CLD, CIRC D (20), HHH,ALFA7, SIGMA, SBETA,XXM, ICPI 
COMMON  IDlIL  , XA,XB,  XC,  I ANG,EP, YC , YR,  llfl TG3, XLBTG3, B I G3, SMAL3 , 333 
COMMON  X3N(5),CCCl,CLt'.eRC,YVY,XM,nERA,SXSin<M,r,X3I00(5),YX3 

-COMMON  P81 7 ,LP,  3ARC.  (-5  ( 3)  rilAllCOCO  I3),LPM,0E,I3R,A3PT,I1I1 

COMMON  BETAN(513),BETAM(5l3j,IJ,LPK,XI!(200),XjJ(200),XaX 
COMMON  AO(7),MPM,MPK,n7ERnD(20),CIRCDI(20),NI3P,NNI3P 
COMMON  I F AG,3XX44 

PAI=3. 141042653  

00  1 1=1,4 

— CC5-I  )=Or 

X3C l)=PAT/t 6. AFLOAT C T-l )*PAI/fl. 

IFTNISP.EO.O)  X3(1)»,0a (PA I/B.+PAI AFLOAT TNNI3P)/64.) 
IFCNISP.CO.O)  X3C2)BPAI/4.-.5ACPAI/4.-2,AXSCn> 

DO  1 .)at,7 

1 CCO-nsCCO-I>AAC.I)ASINCXfiCn  AFLOAT  ( J ) ) 

REtmFw . ..  . ■ ■ - . 


AITKEN  and  DETERM  (both  are  Caltech  subroutines) 
are  not  listed,  although  they  are  included  in  the  deck. 


